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ABSTRACT 


CAUfOftAlu 


An experimental investigation to measure particle size distributions in the 
plume of sub-scale solid rocket motors was conducted. A phase-Doppler particle 
analyzer (PDPA) in conjunction with three-wavelength extinction measurements 
were used in a specially designed particle collection probe in an attempt to 
determine the entire plume particle size distribution. In addition, a laser ensemble 
particle sizer was used for comparative data. The PDPA and Malvern distributions 
agreed in the observed modes near 1 and 4.5 pm diameter (tf). Scanning Electron 
Microscope (SEM) pictures of collected particles were in good agreement with the 
measured Malvern Sauter Mean diameter (djj) of 2.59 pm. Data analysis indicate 
that less than 3% of the total mass of the particles was contained in particles with 
diameter d < 0.5 pm. Therefore, the PDPA, which can typically measure particles 
down to a minimum diameter of 0.5 pm with a dynamic range of 50; 1, 

can be used by itself to determine the particle size distribution. Multiple 
wavelength measurements were found to be very sensitive to inaccuracies in the 
measured transmittances. 




TABLE OF CONTENTS 


I. INTRODUCTION . j 

II. BACKGROUND . 5 

A. PARTICLE SIZE DISTRIBUTIONS . 5 

B. NUMBER, SURFACE/AREA, AND VOLUME/MASS 6 

DISTRIBUTIONS . 

C. SAUTER MEAN DIAMETER . 7 

D. LOG-NORMAL DISTRIBUTION FUNCTION . 8 

E. EXACT SOLUTION FOR THE LOG-NORMAL F(D) .11 

F. MULTIPLE-WAVELENGTH EXTINCTION TECHNIQUE .13 

III. EQUIPMENT DESCRIPTION & SETUP .19 

A. ROCKET MOTOR .19 

1. Propellants . 20 

2. Igniter . 22 

3. Nozzle . 23 

B. PROBE AND TUBES .24 

C. WHITE LIGHT SOURCE AND DIODES .25 

D. PDPA SYSTEM .26 

E. DEFLECTOR .28 

F. DATA ACQUISITION SOFTWARE & HARDWARE .28 

G. MALVERN 2600 . 30 

H. VIDEO CAMERA .30 

IV. EXPERIMENTAL PROCEDURES. 33 

A. PHASE 1 FIRINGS .34 

1. AC-14 10/24/94 . 34 

2. AC-14 11/23/94 . 35 

3. AC-14 12/8/94 . 35 

4. AC-14 1/18/95 . 36 

5. AC-14 2/8/95 . 36 

B. PHASE 2 FIRINGS .37 

1. PS-1 2/23/95 . 37 

2. PS-1 2/25/95 . 37 

3. PS-1 2/27/95 . 38 


































4. PS-1 3/1/95 . 38 

5. PS-1 3/2/95 . 39 

6. PS-1 3/3/95 . 39 

7. PS-1 3/4/95 . 39 

8 . PS-1 3/5/95 . 39 

9. PS-1 3/6/95 . 40 

10. PS-1 3/7/95 . 40 

V. RESULTS & DISCUSSION .43 

A. MOTOR .43 

B. EXPECTED RESULTS .44 

C. MALVERN .46 

D. PDPA .48 

E. THREE-WAVELENGTH EXTINCTION .50 

F. SEM ANALYSIS.53 

VI. CONCLUSIONS .55 

APPENDIX A. MICROPEP OUTPUT.57 

APPENDIX B. PROBE DRAWINGS .61 

APPENDIX C. MIE CODE "MIESCAT3.rOR" .65 

APPENDIX D. MACRO FOR LOTUS l-^3 RELEASE 4 .73 

APPENDIX E. MIE CODE "MIESCATfi-FOR" .75 

APPENDIX F. LINEAR INTERPOLATION IN LOTUS 1-2-3 RELEASE 4 ...83 

APPENDIX G. DATA PRINTOUTS, PLOTS, AND TABLES .85 

LIST OF REFERENCES . 97 

INITIAL DISTRIBUTION LIST. 99 




























LIST OF FIGURES 


1 . Plot of Equations (2-11), (2-12), and (2.13) as PDF#1, #2, and #3 respectively. .. 11 


2. Extinction Coefficients for Three Wavelengths. 16 

3. Average Extinction Coefficients for Three Wavelengths and a Log-Normal 

Distribution. 16 

4. Average Extinction Coefficient Ratios for Three Wavelengths and a 

Log-Normal Distribution. 17 

5. Sub-Scale Motor. 19 

6 - Expected Ideal Pressure-Time Trace for PS-1.24 

7. Spectrum of Oriel Hg-Xe Light Source.26 

8 - Equipment Setup for PDPA and Three-Wavelength Measurements.28 

9- LABTECH Notebook Data Acquisition and Experiment Control Program for 

PDPA and Three-Wavelength Extinction Measurements.29 

10. Malvern Data from PS-1 Firing 3/7/95. 85 

11. PDPA Data from PS-1 Firing 3/2/95.86 

12. PDPA Data from PS-1 Firing 3/3/95.88 

13. PDPA Data from PS-1 Firing 3/4/95.89 

14. LABTECH Notebook Data Acquisition Program Plot of PS-1 Firing 2/25/95. ... 90 

15. Normal Plot of PS-1 Firing 2/25/95 . 91 

16. Bare Light (No Probe) PS-1 Firing 3/5/95.92 

17. PS-1 Firing 3/6/95.93 

18- MATLAB Distribution Plots for "Small Particles".94 

19- Scanning Electron Microscope Pictures of AljOj Collected in the Probe.95 


LIST OF TABLES 


1. Test Propellant AC-14 Specifications. 21 

2. Operational Propellant PS-1 Specifications. 22 

3- Nozzle Selection. 23 

4. Summary of PS-1 Firings Used for Data Reduction. 43 

5. Three-Wavelength Extinction Results. 51 




























I. INTRODUCTION 


Solid rockec motors have used aluminum as a metallic fuel in composite and 
composite-modified double-base propellants for many years. Typically, 5 - 20 % of the 
mass of the propellant is contained in crystalline aluminum particles mixed with 
approximately 70 % ammonium perchlorate (AP) as an oxidizer. The remaining percent 
of the mass is contained in binders, plasticizers, or burning rate catalysts. Metallic fuels 
help performance as measured by an increased specific impulse l^. When metallic fuels 
combust, they react with the oxidizer and form metal oxides. For aluminum, the product 
oxide is A],0,, also called alumina. These compounds are initially in the liquid phase at 
the high temperatures in the motor chamber, but then become solid with the temperature 
drop through the nozzle and into the plume. 

In aluminized propellants the physical process of combustion has been under 
investigation for over thirty years. However, the generally accepted theory is as follows: 
An AP^inder flame starts the process. The 20 - 30 pm aluminum particles melt at 933 K 
on the propellant surface and sinter together forming much larger agglomerates but an 
oxide shell forms around the agglomerate and inhibits further combustion until the 
oxide's melting point is reached at 2327 K. These agglomerates detach the surface and 
enter the flame zone. There the oxide shell melts, allowing inner aluminum to vaporize 
as well as form more oxide. An oxide ’cap" forms on one side. Thus, clouds of burning 
agglomerates, burning unagglomerated droplets, and product oxide droplets exist in the 
motor. When the droplets bum, the vaporized aluminum forms smoke oxide (diameters < 
2 pm) that streams around the agglomerate. The larger droplets with their oxide cap form 
5 - 150 pm residual oxide. Thus, any particle size distribution in the motor is expected to 
bebimodal. [Ref 1-3] 

As the particles pass through the nozzle various processes occur. Collision 
coalescence can cause larger particles to form, whereas, shattering produces just the 
opposite. It has been both measured and modeled that the larger particles entering the 



exhaust nozzle cannot make the rapid turn at the throat, resulting in their being 
concentrated more along the plume centerline. The smaller particles are capable of 
turning with the gas flow and have been observed throughout the plume cross-section. 
Losses occur due to two-phase flow of liquid or solid particles and gas in the nozzle. 
These can be further defined in terms of thermal and velocity lag losses which vary with 
the size of particles. Thus, knowledge of the particle size distribution is important to 
determine these losses. 

Besides determining losses, there are other reasons why particle size information is 
desired. Plume signature and flowfield codes, such as the SIRRM and SPF, require this 
information. Tactically, the AljO, in a missile plume makes the plume very opaque and 
visible, therefore leading to missile detection, tracking, and potentially targeting. For 
example, when viewing a Space Shuttle launch from several miles away, the plume 
emission of the Shuttle Rocket Booster (SRB)'s is highly visible but the huge columns of 
aluminum oxide exhaust are unforgettable. Size distribution information is also 
important for determining the available particle damping in the motor, slag formation, 
heat transfer, and nozzle flow computations [Ref. 1 ]. For example, Blomshield [Ref, 4] 
conducted dry combustion bomb experiments with SRB propellants containing different 
manufacturer's AP to evaluate the variation of the resulting AljO, size distribution 
between them and their impact on slag formation in the SRB's. This work was done as a 
result of the pressure anomaly during STS-54. 

The properties of aluminum oxide are important for use with different panicle 
sizing techniques and calculations. The optical properties of alumina, such as the index 
of refraction, are important for use in optical techniques as well as the previously 
mentioned computer codes. In the plume, both liquid and solid alumina can exist due to 
the range of temperature. The density of AljO, varies with temperature and phase. Reed 
[Ref. 2, p.9] lists the following two formulae for liquid and solid alumina respectively: 


Pl = 2.98-(1.13 X 10-5)(7’-2327) liquid 


( 1 . 1 ) 



Ps = 4.00-(1.16xl0~’)r 


solid 


(1.2) 


where p is density in g/cc and 7 is temperature in K. Salita [Ref. 3, Figure 13] provides a 
plot of the variation of aluminum and alumina densities with temperature as well as their 
ratio. Regarding the complex index of refraction (m = n - ik). Reed [Ref 2] reviewed 
others’ results indicating a range for the AljOj n from 1.65 to 1.76. The complex part k, 
which is the index of absorption, has shown high variability but is almost zero for the 
solid phase in the plume. The liquid phase has a much higher k around 0.5 -1.0 x lO"*. 

The Naval Postgraduate School (NPS) has been conducting research in the area of 
solid rocket motor particle sizing both in the motor and in the plume using various 
techniques for several years. Optical techniques have received the most attention since 
they offer the potential for obtaining panicle information without disturbing the flow, i.e.. 
they are non-intrusive. However, collection probes are sometimes required to limit the 
sample volume because the plume is too opaque for forward scattering measurements. 
Various optical techniques have been applied individually in different regions of the 
motor and plume. Each technique has a limited panicle size measurement range. The 
present investigation used a phase-Doppler particle analyzer (PDPA) in conjunction with 
a three-wavelength extinction measurement in an attempt to determine the entire plume 
particle size distribution. In addition, a laser ensemble particle sizer was used for 
comparative data. The phase-Doppler and three-wavelength extinction measurements 
were made simultaneously in a particle collection probe. A composite propellant with 
20 % aluminum was used in a sub-scale motor to provide the particle-laden plume. 
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II. BACKGROUND 


A. PARTICLE SIZE DISTRIBUTIONS 

Particle size is usually represented by the diameter (in microns) of an equivalent 
size spherical particle. For a collection of particles, size information is often given in the 
form of a size or frequency distribution where frequency is plotted against particle 
diameter. 

There are various types of frequency distributions depending upon the parameter 
of interest. Typically, this is number or mass and the distribution is termed either a 
number or mass distribution. Frequency is another loose term that can represent several 
things. For example, a number distribution could use frequency to mean actual number, 
normalized number, or number percent for a given diameter. 

Due to experimental limitations, frequency is determined empirically for diameter 
intervals where diameter/) is a continuous random variable. Other terms synonymous 
with diameter intervals are classes, class intervals, or size bins. The middle diameter in 
each class is sometimes called the class mark. 

Classification of data leads to representation using a special type of bar chart 
known as a histogram. It uses adjacent rectangles where the base of each rectangle is 
along the abscissa and represents the minimum to maximum diameters of the particles in 
that class. The height of each rectangle corresponds to the ordinate and represents the 
frequency of the particles in that class. For histograms the frequency is usually given per 
diameter. For example, a number distribution would have as frequency the number 
percent in each class divided by the class width. The area of each rectangle then 
represents the number percent in the class. This is a very useful form because a 
probability density function (PDF) p{d) can be fit to the histogram. 

Synonymous terms for probability density ftinction are probability function, 
density function, or distribution function. The PDF can be visualized as the smooth curve 
resulting from connecting the mid-points of the lops of the rectangles as the widths of the 
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rectangles (class widths) approach zero. The area under a PDF equals 1 or 100%, 
meaning the probability P that the diameter is between and 4~ is 100%. Symbolically 
P(-o« < D < +“) = 1, or mathematically 

[2p{D)dD=\ ( 2 . 1 ) 

Since there are no particles with diameters less than zero, the lower limit can be 
set to zero. Similarly, the probability P that the diameter is between two arbitrary sizes a 
and b is just the area under the curve from a to fe or 

P(a<D<fc) =J*p(D)rfD (2.2) 

When the lower limit of the integral above is zero then the integral represents the 
probability P that the diameter is less than b. This has a special name; the cumulative 
probability function or cumulative distribution function/. For a distribution of particles, 
replacing b with d, this function/(<i) represents the percent (or fraction) of particles less 
than diameter d. Mathematically, 

m = -p(D<d) = \'^p{D)dD (2,3) 

B. NUMBER, SURFACE/AREA. AND VOLUME/MASS DISTRIBUTIONS 

Some particle sizing techniques count particles and therefore determine a number 
distribution directly. Others, such as sieving, separate panicles based upon their surface 
area, volume, or weight. Thus it is paramount to know what distribution a given 
technique will yield and to be able to compare these different distributions. Additionally, 
the total number, surface area, volume, mass, or weight of the particles in the distribution 
are used for comparison and analysis. When actual values are of interest as opposed to 
percent (e.g., number instead of number %), then the distribution function often differs 
from the PDF by a factor of the total parameter of interest (e.g., total number of particles 
6 



in the distribution). Take, for example, the number distribution N(d) = N\ PDF where N 
is the total number of particles in the distribution. The PDF is the number percent per 
diameter, and N(,d) is the number per diameter. Integration of N(,d) over all D results in 
N. Integration of N(d) from d, to dj yields the number of particles between d, and d^. 
Since the surface area of a sphere is it (f, the surface or area distribution S(d) 
corresponding to N(d) is S(d) = N(d) xncP where S{d) is the area per diameter. 

Similarly, the volume distribution V(d) corresponding to N(d) is V(d) = N(d) x it d’ / 6 
where V(d) is the volume per diameter. For a constant particle density, the mass 
distribution looks the same as the volume distribution and differs numerically by a factor 
representing the density. Corresponding area and volume or mass PDF's can be found by 
dividing S(d) and V(d) by S and V the total area and volume of particles in the 
distribution, respectively. The constants S and V can be found by integrating S{d) and 
V(d) over all D. Since S(d) and V(d) only differ from their corresponding PDF's by 
constants they have the same shape as the PDF's, just different values. 

C. SAUTER MEAN DIAMETER 

Various diameters are used to compare distributions. They are referred to as 
mean, average, or equivalent diameters. Just as the different distributions are equivalent 
to N(d) times d raised to some power, the mean diameters can be represented for a 
continuous distribution as 

f” N(D)DfdD 

d„--& - (2.4) 

N(P)D‘^dD 


or for a discrete distribution as 

_ 'LN{.d)dfi!d 
- I.N{d)d‘iM 


(2.5) 
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where N{d) is the number per class width 6d [Ref. 5]. One widely used mean diameter is 
the Sauter or Volume-Surface mean diameter dy^- 
D. LOG-NORMAL DISTRIBUTION FUNCTION 

In certain cases, such as ground particles, the histogram has a Gaussian or normal 
distribution shape when the frequency is plotted versus the logarithm of the diameter. 
When plotted on a linear scale the distribution is skewed or stretched to the right. Such a 
distribution is said to be log-normal. 

Like a Gaussian or normal distribution, the log-normal distribution can be 
characterized by two parameters, the first quantifying the central tendency and the second 
determining the spread about this center. For a Gaussian distribution, these two 
parameters are the mean diameter (or arithmetic mean ufn,) and the standard deviation 
c. Similarly, for a log-normal distribution these parameters are the geometric mean 
diameter and the geometric standard deviation Two other diameters are sometimes 
used to describe central tendency: the median and the mode. The median splits the 
distribution in two by number of samples. In probability terms, if all the N sampled 
particles are put in order of increasing size and assigned a number, called the rank (from 
one to AO. then the median diameter is that of the particle with rank (V + l)/2. On the 
other hand the diameter with the highest frequency, i.e., the peak of the distribution, is 
termed the mode diameter. The number of peaks can be specified by referring to the 
distribution as monomodal, bimodal, or trimodal. It is important to note that is also 
the median and the mode diameter for a normal distribution. However, for a log-normal 
distribution d^ is also the median diameter, but is neither the mode nor the mean. For a 
log-normal distribution: 

du = dgexp{-]n^Cg) (2.6) 


dm = £igexp(0.5 In^a*) 


(2.7) 



where dy is the mode diameter and is the mean diameter [Ref. 6]. 

For a log-normaJ distribution of particles, the other representative distributions 
(i.e., number, surface area, or volume/mass) are aJso log-normal with the same <3^. 
However, each distribution has a different d^. To differentiate between them it helps to 
use an additional subscript, For example, d^^ refers to a number distribution whereas d „ 
refers to a mass distribution. In some sources, count with a subscript c is used instead of 
number. The Hatch-Choate Equations relate the various d^'s. One of these is 

=ifj:mexp[-3]n^(o^)] (2,8) 

The following relation is also helpful in relating the various mean diameters to d^. 

=dg»Texpj^^i|—^ln^(aj,)j (2.9) 

= rf*„exp|^^^ln’(o*)j (2,10) 

Since different distributions are of interest to different researchers, rather than discussing 
several d^ 's it is easier to use a mean diameter such as d,j. Note; one aspect of 4,2 can be 
counterintuitive. As one might expect for a given d^„, as <7^ increases the number 
distribution is skewed more to the right and increases. However, if d,^ is held constant 
and Oj increased d^„ decreases and the number distribution shifts left, becoming less 
skewed to the right. The mass geometric mean d^^, however, still increases and the 
volume/mass distribution is more skewed to the right. 

Sometimes log-normal distributions are plotted on log-probability paper since,rtif) 
plots as a straight line with slope o,. The geometric mean d^ is readily identified since it 
corresponds to an ordinate value of 0.5 or 50%. Allthe distributions, i.e., number, area, 
and volume, plot as parallel lines. If the distribution is bimodal rather than monomodal, 
then it has an s-shape between the two lines making up the individual modes. The 
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drawback to the log-probability plot is its unavailability in standard plotting software 
such as spreadsheets. According to the MathWorks however, such a plotting function is 
available for MATLAB In their Statistics Toolbox. 

Several forms of the log-normal distribution function are presented in the 
literature. Three forms are listed below: 


■J^Ogd 


. V2aJ I 


(2.11) 


p(d) = 


y27t log(o^)d 


\0S(d/dg) f 
j2\oe(Cg)_ 


(2.12) 


p{d) = 


ln{<jg)d 


\n{d/dg) 
J2 ln(Oj,). 


(2.13) 


Some sources use a version of one of the three equations above without the d in the 
denominator. [Ref. 6-10] 

It is therefore important to determine the correct form. Three criteria are used to 
check them: 

1. Integration of the PDF (area under the curve) from d = 0 to -t-®® must equal one 
(i.e.,/(+»°)= 1). 

2. Using^ti), the ratio of d^^^ (d whert fid) = 0.84 or 84%) over d^^ {d where 
fid) = 0.50 or 50%) must equal c^. 

3. The function must behave as depicted in references for different values of d^ 
and Oj. 

MATLAB can be used for analysis. The functions QUAD and QUADS performed 
numerical integration, although an exact solution fotfid) is available for Equations (2.11) 
and (2.13) as described in the next section. Figure 1 shows a plot of Equations (2.11) to 
(2.13). All three equations without the d in the denominator fail criterion 1. However, it 
should be noted that some sources integrate the PDF with respect to ln(<f) vice d. This is 


10 



the same as integrating a PDF, with a t/ in the denominator, with respect to d since 
<i[ln(D)] = {\ID)dD in differential form. As is, Equation (2.12) fails criterion 1 as is 
readily apparent by inspecting the figure. Equation (2.11) fails criterion 2. Equations 
(2.11) and (2.12) fail criterion 3. Equation (2.13) passes all criteria and thus was used in 
the analytical work. Note however that Equation (2.11) is a PDF and fits in as Kerker 
[Ref. 6, p. 357] describes "a family of logarithmically skewed distributions”. [Ref. 6] 


Possible Log-Normal PDFs 

Dg=3. SIGMAg=1.5 
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Figure 1. Plot of Equations (2.11), (2.12), and (2.13) as PDF #1. #2, and #3 
respectively. 


E. EXACT SOLUnON FOR THE LOG-NORMAL F(D) 


The general form of the log-normal distribution function is 































where p and a are constants [Ref. lOJ. Most often, for solid rocket motor particle sizing 
P = infop and a = ln(dp as in Equation (2.13). Although numerical integration 
techniques can be used to obtain an approximate^<f), there is an exact solution using the 
error function erf(jt). 

erf(j<) = -L f“exp(-H)* (2.15) 

Jn ■’0 


and using Equations (2.3) and (2.14) 

w-J' 


0 V27t PZ) 


ln(D) - a 

. V2P . 


By letting 


ln(P)-a 

VIp 




(2.16) 


(2.17) 


di = 


yipD 


dD 


and the limits of integration are 


lim ln(P)-P ln(<f)-a 

VIp J” VIp 


upper limit 


(2.18) 


(2.19) 


ln(P)-tt 

. V2|3 


lower limit 


( 2 . 20 ) 


Therefore, 
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(2.21) 






VifJ 


J“exp(-/=')dr= Jn 


(2.23) 


and exp(-/^) is an even function, i.e., symmetric about the ordinate, then 

^Jlexp(-a.<>=^Jrcxp(-a^.4 


Using this and the definition of the error function, then 


M = r + ierf 


ln(rf)-a 

VIp . 



(2.24) 


(2.25) 


[Ref. 11]. Thus for solid rocket motor particle sizing 



+ erf| 


]i\(,d/dg) 1 ] 

V2 ln((j^)Jj 


(2.26) 


F. MULTIPLE-WAVELENGTH EXTINCTION TECHNIQUE 

Cashdollar [Ref. 12] discusses the use of light extinction measurements for 
extracting particle size distribution parameters. Light entering a cloud of particles is 
absorbed, scattered, and transmitted. Transmittance T is the ratio of transmitted intensity 
to incident intensity. Light extinction is the sum of absorption and scattering. When Mie 
theory is applied to light scattered by single spherical particles the extinction coefficient 
Q can be determined as a function of wavelength X., particle diameter d, and complex 
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index of refraction m~n-ik(k is the absorption index). Bouguer’s law provides the 
transmittance through a nionodisperse (single size) distribution of particles; 


T = exp I 


'iQC„L \ 

, J. 


(2.27) 


T - transmittance 

Q = dimensionless extinction coefficient 

C„ = mass concentration of particles (mass of particles per volume of gas) 
L = path length through the cloud of particles 
p = particle density 
d = particle diameter 


The volume concentration of particles (volume of particles per volume of gas) can be 
related to C„ by 

C..^ (2.28) 


To obtain accurate values for extinction, it is necessary to have the detector not see 
forward scattered light. An approximation for the required detector viewing half angle 
to satisfy this requirement is [Ref. 12] 

01 <7.0^ (2.29) 


Dobbins revised the Bouguer transmission law to apply to a polydispierse system (many 
sizes) [Ref. 12]: 



where Q is the average extinction coefficient, found for a continuous distribution by 


14 



_ r Q(D)N{D)D^dD 

Q = - 

\ N{D)D^dD 


(2.31) 


or for a discrete distribution by 

- _ 'LQ{d)N(d)d^&d 
^ l.N(d)d^M 


(2.32) 


Since Q depends upon the distribution, then for a log-normal distribution it is a function 
of <fj2. Oj, X, and m. Taking the natural logarithm of both sides of Equation (2.30) and 
ratioing two equations for different wavelengths results in 


ln7-(Xi) g(X|) 

ln7'(X2) Q{X2) 


(2.33) 


This is the basis of the multiple-wavelength technique. Transmittances for different 
wavelengths of light passing through the same volume of particles over identical path 
lengths are experimentally measured. A Mie code is used to generate Q's for the set 
wavelengths and for various values of tfjj. Ug, and m. Then, using Equation (2.33) the 
ratios can be compared and further iterations performed as necessary to determine the 
best values of d^^, and m. The more wavelengths used the more ratios to compare and, 

presumably, the more accurate the method. The number of ratios for a set of n 
wavelengths is equal to the sum of integers less than n (e.g., for six wavelengths there are 
1+2 + 3 + 4 + 5*15 ratios). Note that the number of independent ratios is only n - 1 
(e.g., five for six wavelengths). The accuracy of the method is also increased by 
spreading the wavelengths across the spectrum [Ref. 13). To determine dj^, and 
m = n - ik (4 variables) a minimum of four independent ratios are required. Figures 2, 3, 
and 4 show plots of the Mie code Q's. Q's. and Q ratios for the three wavelengths used in 
the experiments and an assumed log-normal distribution. Notice from Figure 4 that as djj 
increases the Q ratios flatten out. For this reason the technique does not work for 
distributions with lots of mass in large particles. 
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Figure 2 Extinction Coefficients for Three Wavelengths. 



Normal Distribution. 
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Figure 4. Average Extinction Coefficient Ratios for Three Wavelengths and 
Log-Normal distribution. 


























III. EQUIPMENT DESCRIPTION & SETUP 


A. ROCKET MOTOR 

For all of the test firings, a sub-scaJe motor was used as shown in Figure 5. The 
motor consists of a stainless steel cylindrical chamber with a 2" inside diameter, 0.625" 
wall thickness, and overall length of 9.25 '. The head of the chamber is closed off by a 1" 
thick end plate with rubber o-ring seal and secured by six bolts. The nozzle end of the 
chamber has a slightly larger inside diameter to allow the nozzle insert to slide in place, 
and is sealed by another end plate and six bolts. A simple converging-diverging nozzle 
insert is used. Gas leakage around the nozzle is prevented by the use of two o-rings fit 
into grooves around the exterior of the nozzle insert. A pressure transducer and pressure 
relief fitting are screwed into the chamber. Prior to motor firing, a hydraulic dead-weight 
tester is used to determine the calibration slope and zero for the pressure transducer. The 
pressure relief fitting consists of a steel elbow with a 1000 psia burst disk. The propellant 
is bonded to the head-end plate and inner chamber wall using Permatex High Temp RTV. 
The igniter consists of a 3/8" steel bolt, described below, that is screwed into the optimum 



Figure 5. Sub-Scale Motor. 
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of four available holes angled toward the head of the chamber. The remaining holes are 
filled with standard bolts. All o-rings are lubricated with Dow High Vacuum Grease. 

1. Propellants 

Two different propellants provided by the Air Force Phillips Laboratory were 
used. A 1" thick and approximately 2" diameter propellant grain was cut from the stock 
using a large knife, then the desired cutouts were made using a three-piece steel cutting 
tool and a press. The tool could make 0.6" I.D., 1.4" O.D., 1.5" I.D. and 2" O.D. cutouts. 
The grain geometry was chosen to optimize the burning surface area for a neutral bum 
with a chamber pressure of at least 250 psia lasting for about one second. RTV was used 
to bond each grain in the motor as well as to inhibit certain surfaces from burning. 

For testing of the experimental setup, especially the three-wavelength apparatus, 
several firings were conducted using a propellant (AC-14) that was on hand. Its 
composition is shown in Table 1. It was undesirable for data analysis due to the silicon 
which caused oxides of silicon or aluminum together with mullite in the exhaust. To 
achieve the desired performance several grain configurations were tried. Initially, a single 
1" thick end-burning grain (2" O.D.) was tried, then two 1" thick grains were tried, one 
end-buming (1.4" O.D.), the other radial- and end-buming (1.5" I.D., 2" O.D.). Finally, 
good performance was obtained using two 1" thick radial- and end- burning grains. The 
first had a 2" O.D. and a 1.5" I.D. (0.25" web). RTV was used to inhibit the downstream 
end and bond the grain to the chamber wall just downstream of the igniter hole. The 
second grain was cut from the center of the first. It had a 1.4" O.D. and a 0.6" I.D. (0.4" 
web). RTV was used to bond the upstream surface to the endplate and to inhibit the outer 
radial surface. 
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ID 

AC-14 

Composition: AP 

67.15% 

A1 

12.0% 

R45M 

10.14% 

Si 

6.0% 

Dioctyl Adipate 

3.91% 

EPDI 

0.78% 

Triphenyl Bismuth 

0.02% 

Density (Ibm/in’) 

p„ = 0.06 

Burning rate pressure coefficient 

a = 0.0086 

Burning rate pressure exponent 

n = 0.566 

Characteristic exhaust velocity (ft/s) 

c' = 4957 


Table 1. Test Propellant AC-14 Specifications. 


As previously mentioned aluminized propellants have been the subject of particle 
sizing for the past thirty years. Therefore, it was desired to obtain such a propellant since 
various parameters are better known (e.g., particle index of refraction) and this enables 
comparison to past research. The USAF Phillips Lab provided four propellants, all of 
which were class 1.3, conventional composites. The one selected for this thesis research 
was PS-1. Its composition (by mass %) and characteristics are listed in Table 2. Initially, 
a single 1" thick end-buming grain (2" O.D.) was tried; however, the required small 
nozzle throat diameter and propellant's high aluminum content resulted in clogging of the 
throat. A single perforated end- and radial-burning grain was then chosen and used for 
the remainder of the research. It was 1" thick with a 2" O.D. and a 0.625" I.D. 
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ID 

PSI 1 

Composition: AP (coarse) 

50% 200nm 

AP (fine) 

20% 

A1 

20% 

HTPB 

10% 

Density (Ibm/in’) 

p„ = 0.0667 

Burning rate pressure coefficient 

a = 0.0328 

Burning rate pressure exponent 

n = 0.376 

Characteristic exhaust velocity (ft/s) 

c' = 5119 


Table 2. Operational Propellant PS-1 Specifications. 


2. Igniter 

Each pyrotechnic squib igniter was processed as follows: First, a hollow 3/8” 
steel bolt had two holes drilled in the bolt head to accomodate the lead wires. Two 
approximately 8" long copper lead wires were then cut and some insulation was stripped 
off both ends with an X-Acto knife. The wires were passed through the holes and then a 
small piece of nickel-chromium wire was soldered (using ruby flux and pure solder) to 
the ends leaving about a 1/16" gap. Next, the wire bridge was pulled inside the cavity in 
the bolt and the bolt head was epoxied on the outside to secure the wires. After the epxjxy 
had dried the cavity was filled with BKNO, and a paper hole from a hole punch was 
glued on the end with Elmer's Glue. Continuity checks were performed at various times 
in the process to verify that the wires were not shorted by touching the inside of the bolt. 

For firing the igniter, lead wires were connected to a 12 volt battery with battery 
charger and a firing button in the control room. When the firing circuit was completed by 
depressing the button, the heat dissipated in the nickel-chromium bridge wire flashed the 
BKNO,. The burning matter burst through the paper hole and impacted the grain, thereby 
initiating propellant combustion. 
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3. Nozzle 


A copper nozzle insert was selected to obtain reasonable chamber pressure. The 
following steady-state mass conservation equation was used: 




~ chamber pressure 
Aj = burning surface area 
A, = nozzle throat area 
Pj, = propellant density 
a = burning rate pressure coefficient 
c = characteristic exhaust velocity 
r|,. = characteristic exhaust velocity efficiency 
= gravitational constant (for English units) 
n = burning rate pressure exponent 


(3-1) 


The MICROPEP computer program [Ref. 14] was run with the propellant ingredients to 
obtain and c. See the output in Appendix A. Based upon previous results, the 
efficiency was assumed to be 0.9. Then, the available nozzle throat areas were 
substituted to calculate p^. Due to clogging of the throat from the high aluminum content, 
the largest A, that would yield p^ > 250 psia was selected. The copper nozzles were 
converging-diverging and Table 3 shows the selection that was available. 

For initial testing with propellant AC-14 several nozzles were used. The 


d,(ii\) 

0.207 

0.340 

0.239 

0.165 

0.452 

0.216 

0.270 

0.200' 

0.208 

0.268 

0.286 

d, (in) 

0.378 

0.440 

0.365 

0.230 

0.662 

0.395 

0.808 

0.368 

0.425 

B1RB1 


e = AJA, 

3.335 

1.675 

2.332 

1.943 

2.145 

3.344 

8.956 

3.386 

4.175 

3.481 

1.268 

Con/Di 

(deg) 

45/15 

45/15 

45/15 

45/15 

45/15 

45/15 

??? 

??? 

??? 

45/15 

??? 


Table 3. Nozzle Selection Where d, is the Nozzle Throat Diameter, d, is the Nozzle 
Exit Diameter, £ is the Nozzle Area Ratio, and Con/Di are the Nozzle Converging and 
Diverging Half-Angles Respectively. 
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d, = 0.165" nozzle was tried with the single end-burning grain and resulted in twice the 
expected p^, probably due to clogging. Next, the first two-grain configuration was paired 
with the d, = 0.239" nozzle and resulted in a progressive initial burn and a long low 
tail-off. Finally, the second two-grain configuration was paired with the d, = 0.268" 
nozzle with good results. 

For the PS-1 firings, two nozzles were used. The d, = 0.207" nozzle was tried 
with the single end-burning grain but clogging caused a highly progressive bum. For use 
with the single end- and radial-burning grain, the d, = 0.239" nozzle was chosen and 
Figure 6 shows the expected ideal performance. The actual performance resulted in a 
higher chamber pressure of about 350 psia. Post-firing inspection of the nozzle revealed 
a fairly uniform layer of AljO, reducing d, to about 0.15". 



Figure 6. Expected Ideal Pressure-Time Trace for PS-1. 


B. PROBE AND TUBES 

A particle collection probe that consists of a narrow, hollow steel box with a 
conical tip for swallowing a portion of the plume exhaust was previously developed [Ref. 
20]. The design of the probe was to swallow the normal shock, forming weaker, internal 
oblique shocks that would not break up the particles. It was positioned to avoid Mach 
disk locations. Windows on either side had been used with optical sizing techniques and 
a filter on the downstream end had been used to collect particles for scanning electron 
microscope (SEM) analysis. For the current uses, the windows were removed and steel 
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lubes added on the sides to enclose the PDPA beams and to vertically enclose the Hg-Xe 
beam described below. The insides of the tubes were spray-painted flat black. To protect 
from bumthrough, an additional steel plate was welded on the funnel side tube. To 
ensure that the PDPA lenses would not be coated by exhaust, circular plate glass was 
taped in the ends of the side tubes. A nitrogen line with a sonic choke was connected to 
the probe and provided an ejector flow around the probe tip where it exhausted into the 
probe body. The back pressure in the probe was regulated to provide and approximate 
isokinetic sampling condition. For some runs particles were wiped from the inside of the 
box for SEM analysis. RTV was put on parts of the front of the probe body for thermal 
protection. See the diagrams in Appendix B. 

C. WHITE LIGHT SOURCE AND DIODES 

An Oriel Model 66002 Hg-Xe arc lamp with a f/1 condensing lens and 200 W 
power supply was used as the light source. It was mounted below the lab table and a 
magnetically held mirror was positioned to reflect the beam up through a 1/2" diameter 
hole in the table. From the top of the table, tubing enclosed the beam up into the probe 
and above the probe to the diode box. Both tubes were spray-painted flat black inside. 
The lower tube had a 0.5" O.D. and 7/16" I.D. and a collar adapted it to the probe's lower 
vertical tube O.D. A larger upper tube (1" O.D.. 7/8" I.D.) was used to cut down on 
internal reflection and was fit over the probe’s upper vertical tube. The diode box 
contained an adjustable pinhole, two beamsplitters, three EG&G Model HUV-4000B 
diodes, each with a Melles Griot narrow band interference filter (0.01 pm bandwidth) to 
measure wavelengths 0.3130.0.5461. or 0.8500 pm. neutral density filters, and an 
NPS-built preamplifier for the diodes. The pinhole and neutral density filters were 
adjusted to get as high an output voltage as possible but less than 10 volts (LABTECH 
Notebook limitation; the diodes saturated near 15 volts). The total length from the 
condensing lens to each diode was about 80". From the probe center to the diode pinhole 
was about 36". Therefore, since the pihole was open to about </ = 0.19" the detector 
viewing half-angle was about 0. IS", easily satisfying the requirement of Equation (2.29). 
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For one run a light chopper was used between the arc lamp and mirror Figure 7 (lower 
curve) shows the spectral irradiance from 0.200 to 0.900 pm of the light source and the 
bandpass filter locations. 

D. PDPA SYSTEM 

The Aerometrics Phase Doppler Particle Analyzer (PDPA) system consists of a 
laser transmitter, receiver, doppler signal analyzer (DSA), and computer digital signal 
processor (DSP) and software. A Lexel Model 95-4 argon ion laser, tuned to 0.5145 pm 
and operating at 2 W, with beam waist adjustor, polarization rotator, and steering mirrors 
provided the input beam to the transmitter. The transmitter uses a Bragg cell to split the 
incoming laser beam into two beams of equal intensity, a Zero order beam and a 40 MHz 
shifted First order beam, then bends the beams with prisms and mirrors and converges 
them with the exit lens. The beams converge until they cross, forming an ellipsoidal 


Gaussian volume, approximately 10" in front of the transmitter. When a particle passes 
through the crossing volume, light is scattered. For panicles with t/> 0.5 pm this is 



Figure? Spectrum of Oriel Hg-Xe light source [Ref 15,p. 1-41,fig 8] 
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scattering in the Mie and geometric optics regimes of scattering. To simplify processing 
the system assumes only the geometric optics regime applies. The receiver, consisting 
essentially of optics and photo multiplier tubes, is positioned at a predetermined optimum 
angle to collect a portion of the scattered light. Aerometrics ran simulations, given 
estimated particle indices of refraction for different motor locations, to determine this 
optimum angle for a linear phase response over the range of particle diameters. They 
found that a scattering angle of 130° (backscatter or reflection) was best for 
through-the-chamber particle sizing and a scattering angle of 45° (forward scattering or 
refraction) was best for through-the-plume measurements as presented in this paper. The 
receiver was mounted on a support at a 45° angle and approximately 8” from the crossing 
volume. The horizontal tubes of the probe were designed to fit around the collars of the 
transmitter and receiver. [Ref. 16] 

Data acquisition and processing was done with a personal computer (PC) and the 
DSA. The DSP was located in the PC and was used to quickly perform discrete and fast 
Fourier transforms on the data. The DSA performed several functions. It provided the 
input to the transmitter to create the two beams and conducted calibration and data 
processing. The software created the user interface through menus for setup, calibration, 
data acquisition, processing, display, and storage. The PDPA was externally triggered to 
take data through the LABTECH Notebook program run on a separate PC as described 
below. [Ref. 16] 

The PDPA software rejects particles for various reasons. The most important 
ones for this paper are multiple panicles and undersize particles. Only single panicles 
can be accurately accounted for, therefore if the system detects multiple particles that 
sample is rejected. Panicles with diametersas small as approximately 0.2 pm are detected 
but if less than 0.5 pm they are rejected by the system. No information is available for 
panicles less than 0.2 pm since the scattered light intensity is loo low to be detected. The 
PDPA has Probe Volume Correction (PVC) that allows for correcting the measured 
distribution for the smaller panicles. However, at least 100 valid samples must betaken 
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Figure 8 shows the experimental setup for the PDPA and three- wavelength extinction 
measurements in the probe. 



Figure 8. Equipment Setup for PDPA and Three-Wavelength Measurements (Deflector, 
Video Camera, and Supporting Hardware Not Shown). 


E. DEFLECTOR 

A pneumatically actuated, computer controlled steel deflector plate was positioned 
between the motor nozzle and probe. Its purpose was to protect the probe and only allow 
exhaust entry into the probe during steady burning (i.e., the plateau of the chamber 
pressure-time trace for a neutral burning grain). Some firings burned holes through the 
plate, requiring repair by welding and reinforcement. 

F. DATA ACQElSrnON SOFTWARE & HARDWARE 

A 486DX/33 MHz PC controlled the experiment through a Windows''^'^ 
application known as LABTECH Notebook. Icon blocks, representing analog and digital 
input or output devices, with connecting arrows were used in the software to visualize the 
control of equipment and the flow of data. See Figure 9 Each icon block could be 
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opened to set various parameters. When the data acquisition file was nin the following 
events were executed: 

1 The diode and pressure transducer outputs were recorded for the entire run (20 
seconds), 

2 The nitrogen ejector flow was turned on from 2 to 18 seconds, 

3. An audio signal was sent to the videorecorder at 0.5 seconds as a time 
reference, 

4. When the motor chamber pressure reached 100 psia a 1 second delay occurred, 
then the deflector was triggered to slide out of the way (for 0 7 seconds) and 
the PDPA was triggered to take data 0.25 seconds later, and 

5. At the termination of the run the clock, chamber pressure, three diodes, and 
deflector position data were output to a separate file. 



Figure 9. LABTECH Notebook Data Acquisition and Experiment Control Program 
for PDPA and Three-Wavelength Extinction Measurements. 
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Keithley MetraByte hardware provided the interfaces between the computer and 
the equipment. The 24-Bit Parallel Digital I/O Board (PIO-24) and the Analog & Digital 
I/O Board (DAS-16F) were internal to the computer. Externally, the PIO-24 was 
connected to a 24 Channel Relay Output Board (ERB-24), The DAS-16F connected to 
the Screw Terminal Accessory Board (STA-16). 

G. MALVERN 2600 

Ensemble particle-sizing was done using the Malvern 2600. This device consists 
of a laser transmitter and receiver and computer for control and data acquisition. A 9 mm 
diameter He-Ne laser beam with wavelength X = 0.633 jim exits the transmitter and is 
scattered (Fraunhofer diffracted) by the particles based on their size. The scattered light 
then passes through a Fourier transform len and is focused on 32 semicircular detector 
rings to measure the intensity of light at different solid angles up to a scattering angle 0 of 
14°. Due to the physical limitations of positioning the Malvern, it was not possible to 
simultaneously conduct Malvern, PDPA, and three-wavelength extinction measurements. 
Therefore, the Malvern was used to take a sample through the plume at the same 
longitudinal position 6" downstream of the nozzle as the main PDPA measurements. It 
was triggered using a separate LABTECH Notebook program to take 50 sweeps of data in 
0.4 seconds, starting 1 second after the chamber pressure reached 100 psia. [Ref. 17] 

H. VIDEO CAMERA 

A video camera was positioned above the setup and framed on the plume region. 
The signal was sent live into the control room where it was monitored and recorded on a 
Sony S-VHS VCR. Text annotation and a time reference were displayed on the screen, 
however the time signal was not synched with the LABTECH Notebook clock. 

Therefore, an approximate time hack was sent from the data acquisition program in the 
form of an audio signal that would spike the VU meter on the VCR. The display was 
used to measure the plume dimensions, determine probe tip position relative to the plume 
centerline, and safely observe the progress of the run. Although 1" separated holes in the 


30 





lab table were typically visible, parallax prevented their direct use in measuring motor 
nozzle to probe tip posiiion. 
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IV. EXPERIMENTAL PROCEDURES 


Preparation for each firing required several steps as discussed below. Since the 
igniters required extra time to fabricate, they were normally made in advance during 
non-firing days. Each propellant grain was cut, bonded to the motor case and allowed to 
cure overnight. On the day of the firing the probe and side tubes were assembled and 
coarsely aligned on their mount using pencil marks of the final previous position as a 
guide. The argon laser, with its water cooler and fan, and the PDPA were then energized 
to conduct fine alignment. Cigar smoke, deodorant spray, or humidifier water vapor were 
used in the probe to see the beam and test the PDPA. The white light source was also 
turned on and its beam checked using a white card at various spots. The crucial check 
was a clean, collimated beam exiting the top vertical tube of the probe. Next, the diode 
box was coarsely aligned and the extension tubes were put in place. The diode box was 
opened and the filters and beam splitters cleaned and secured. The incoming beam 
exiting the pinhole was positioned on the center of each diode. The cover was replaced 
and the diode outputs verified. The propellant surface was scraped to provide a fresh 
surface to aid in ignition. The motor was assembled without the igniter and mounted at 
the proper longitudinal, lateral, and vertical position in front of the probe tip. The video 
camera was turned on and adjusted. The probe Nj ejector flow valve was activated to 
purge the line and then connected to the probe. The deflector was aligned and the 
activator air was turned on. All computers were booted and data acquisition software 
loaded. The videotape recorder, monitor, and screen annotator were turned on and set up. 
The white light beam was covered temporarily so "no light" diode readings could be 
taken. Several dry runs with the LABTECH Notebook program were performed and a 
manually actuated voltage to simulate the pressure transducer was input to fully check 
system performance and data acquisition. Old files were deleted to ensure adequate hard 
drive space for the 0.5 MB data file. Finally, protective shields and plastic covers secured 
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with duct tape were attached to protect the PDPA from exhaust products. The igniter was 
installed and battery/charger hooked up. 

The firing sequence started by clearing the perimeter area of personnel. Next, the 
warning horn was sounded, VCR set to record, PDPA set to acquire data when externally 
triggered, and the LABTECH Notebook program started. The final event was a short 
countdown ''3-2-1-FIRE" and the ignition button pressed and held until the firing was 
initiated or the program was terminated in the event of a "no-fire". 

Initial post-firing procedures consisted of stopping PDPA data acquisition and 
saving the data, stopping the VCR, and securing the firing circuit. The laser, arc lamp, 
diodes, and video camera were turned off. The battery charger was unplugged and N^/air 
valves closed for the ejector and deflector. If SEM samples were to be taken, then they 
were obtained by wiping from the inside of the back of the probe with a surgical gloved 
finger dipped in acetone. Lenses, mirrors, and other glass or optics were wiped with 
methanol or spray-cleaned with compressed air. All exposed areas were wiped with 
acetone to try to prevent corrosion from the exhaust. The probe and all tubes as well as 
the motor were disassembled and cleaned with water or acetone. 

A. PHASE 1 FIRINGS 

Prior to the operational propellant PS-1 firings, several test firings were conducted 
with AC-14 to check the setup and to make any necessary modifications. Measurements 
using the PDPA and original probe were made and drawings made of the required 
changes and additional hardware. Several months were required to get the probe 
modified (new hole cut and vertical tubes added) and the PDPA side tubes machined. A 
collimation tube with pinhole and lens holders was manufactured for use with the arc 
lamp. Additional probe tips (spares) were also made in case some were damaged during 
firings. 

1. AC-14 10/24/94 

This first firing was to check whether an end-burning grain and the £/, = 0.165" 
nozzle would yield the desired 250 psia without clogging. It also was an opportunity to 


34 



test the PDPA through the newly modified probe with side tubes. The probe tip was 
positioned 5.5" (24 d,) from the nozzle and 1” (4 rf,) off of centerline. Unfortunately, the 
nozzle clogged with molten aluminum and silicon forcing the chamber pressure up to 600 
psia. The data acquisition program, N2 ejector, deflector, and PDPA appeared to work as 
desired, although the trigger times needed adjustment. 

2. AC-14 11/23/94 

The second firing was an opportunity to try a different grain configuration and 
nozzle combination. Two grains, one end- and the other end- and radial-burning 
(described in Chapter HI), were used with the d, ® 0.239" nozzle to ideally achieve 350 
psia for about 1 second. However, a somewhat progressive bum resulted in a chamber 
pressure of 470 psia. Clogging was still a problem. The probe was positioned 5.5" (15 
d,) downstream and 1.5" (4 d^) laterally. The PDPA did not acquire any meaningful data. 

3. AC-14 12/8/94 

This was the first run with the three-wavelength extinction apparatus in addition 
to the PDPA. A small diameter tube connected the top vertical tube of the probe to the 
diode box. The arc lamp was positioned on a stand next to the lab table. The collimation 
tube was attached to the condenser housing on one end and on the other end to a box 
containing a 45° mirror beneath the lower vertical tube of the probe. It was discovered 
that the neutral density filters and beamsplitters were blocking UV light (0.365 or 0.313 
pm), so the UV filter was moved to the lowest diode in the box and the neutral density 
filter at the box entrance was removed. This also required changing neutral density filters 
in front of the other diodes until suitable output levels were attained. The other two-grain 
configuration for AC-14 described in Chapter III was used with the d, = 0.268" nozzle to 
ideally achieve 300 psia for 1 second. The actual pressure peaked at 420 psia and was 
regressive burning for approximately one second. Since this grain and nozzle 
combination worked quite well it was used for all subsequent runs with AC-14. The 
diode outputs initially started to drop as expected, then saturated, then decayed to a new 
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full light value. This was interpreted as a vibration or shock problem resulting from 
having the lamp, collimator, and mirror too close to the motor. The Nj ejector line had 
been damaged on previous mns and was burned off this time. 

4. AC-14 1/18/95 

For the fourth test firing several changes were made. A hole was drilled in the 
table and the white light source was moved below the table. An adjustable minor on a 
post with a magnetic mount was added below the table and replaced the minor in the box. 
It was also determined that a more intense and clean beam resulted without the secondary 
collimator, so it was no longer used. The light was collimated instead simply by the lamp 
housing lens and the 80" distance from the lamp to the diodes. The probe had a new Nj 
ejector connector hole drilled on the opposite side for better protection. The pressure¬ 
time trace contained a 1 second neutral bum at 400 psia and the diodes responded 
properly although at a questionable rate. It was noted that the diode outputs dropped 
during Nj ejector activation, probably due to circulation of stray particles. The vibration 
problem appeared to be solved. A hole was burned in the side fuimel which connected 
the probe body to the PDPA detector. 

5. AC-14 2/8/95 

The side funnel was repaired and had a protective plate added in front. The diode 
response question was handled in several ways. Unnecessary amplifiers that had a 1 Hz 
filter were disconnected. The diode outputs were then increased by adjusting trim pots 
and swapping out neutral density filters. The sampling rate for the diodes was increased 
to 300 Hz (this also applied to the pressure transducer and deflector trigger). Finally, a 
light chopper was used to produce a 15 Hz square wave and the diodes had no difficulty 
measuring the alternating light conditions. It was also considered that a sluggish deflector 
could have caused the appearance of slow diode response. Therefore, the video replay 
was used to verify that the deflector was operating as commanded as accurately as could 
be measured. This was the final AC-14 firing because everything seemed to work 
properly. The acquisition run time was decreased from 30 to 20 seconds. The deflector 
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retraccion interval was increased from 0.65 seconds to 1 second. The probe tip was 
located 5.5" (11 d,) from the nozzle and 1" (5d^) off centerline. The Nj ejector flow was 
regulated at 300 psia, but may have been too high (non-isokinetic flow) as evidenced by 
the diode outputs dropping prior to the deflector retracting. In addition, good 
transmittance data was obtained although the UV diode nearly reached T = 0. This was 
probably due to too large a regulated Nj pressure. The pressure-time trace showed an 
average of 350 psia for about 1 second. Good PDPA data was obtained (97 valid 
particles with a d^^ = 10.3 ftm). 

B. PHASE 2 FIRINGS 

Since the AC-14 firings had enabled fixing most of the experimental problems, 
the operational firings with PS-1 were able to proceed at a much faster pace. The new 
propellant had good burning characteristics, ignited easily, and did not produce the erratic 
"puffing" of exhaust that was typically observed with other propellants. It also was easy 
to cut and scrape. However, the high aluminum content of the propellant produced a very 
difficult to remove AljO, layer in the nozzle. 

1. PS-12/23/95 

As with AC-14 an end-burning grain was initially tried for simplicity and 
assurance of neutral burning. The d, = 0.207" nozzle was used and ideally should have 
resulted in a chamber pressure of 224 psia for 4 seconds. However, nozzle clogging was 
a major problem and a highly progressive burn resulted with a peak chamber pressure of 
430 psia. The probe, PDPA, and three-wavelength apparatus were not used since this 
first run was for plume visualization and measurement. The clogging problem caused the 
plume to be irregular, so this run was not as useftil as desired. 

2. PS-l 2/25/95 

A new calibration curve was generated for the pressure transducer prior to this 
run. A single perforated grain and the larger d, = 0.239" nozzle were selected and used 
for the rest of the work. Ideally an approximately neutral bum near 300 psia for 1.5 
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seconds would occur. The deflector and PDPA triggering were changed to the intervals 
0.25 - 1.25 and 0.45 - 1.25 seconds after reaching 100 psia respectively. The Nj ejector 
regulated pressure was decreased to 100 psia. The motor was positioned for a 
nozzle-to-tip distance of 6" (16 rf,) and a lateral offset of 1.2” (3 d^). The diode box 
pinhole was opened slightly creating a beam with a diameter of about 0.29". The narrow 
top tube from the probe to the diode box was replaced with a much larger 1" O.D. tube to 
cut down on internal reflection, improving the beam. The actual chamber pressure 
peaked at 400 psia. The PDPA did not measure any particles. The transmittances looked 
good, however they dropped to a very low value. Also, because the transmittances started 
dropping early it was felt that the Nj pressure was again too high. From the video replay 
it appeared that the tip was located in the smoke periphery of the plume as desired. The 
plume width was estimated to be 3" with the highly visible portion about 1.5" wide. 

3. PS-12/27/95 

This run was to be on centerline. Since only large particles {d > 0.5 nm) were 
anticipated, only the PDPA was used. Concern for being too close and melting the tip 
drove the nozzle-to-tip distance up to 12" (33 </,). To reduce the tip's exposure time to the 
hot plume the deflector trigger interval was shortened to 0.25 - 0.95 seconds after ICX) 
psia chamber pressure was reached. The actual run had an essentially neutral bum at a 
chamber pressure of about 3(X) psia for 2 seconds. Unfortunately, the PDPA only 
measured 8 valid panicles with a d^^ = 12.0. It was decided that the three-wavelength 
measurements should always be taken to aid in evaluating what the PDPA measured. The 
nozzle throat diameter was reduced by the AljO,, and was measured after the run to be 
0.15". This was typical of all the subsequent runs. 

4. PS-13/1/95 

From the two previous runs it appeared that the deflector and PDPA were being 
triggered early, so the intervals were again changed to 1 - 1.7 and 1.25 - 1.7 seconds after 
100 psia respectively. The nozzle-to-tip distance was significantly reduced to 4" (11 d^ 
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and this run was on centerline. No PDPA data was collected, possibly due to a computer 
problem after an automatic laser shutdown occurred after an overheat condition with the 
water cooler. The probe tip melted, probably 0.3 seconds after deflector retraction- 
particles were collected for SEM analysis. 

5. PS-13/2/95 

The nozzle-to-tip distance was set to 6" (16«/,) and 0.2” (0.5 laterally. The 
chamber pressure peaked near 380 psia had an essentially neutral bum. Transmittance 
data from this and the previous run showed that the ejector flow was still too high. Good 
PDPA data was acquired with 44 valid particles. Particles were collected again for SEM 
analysis. The tip melted, probably 0.1 seconds after deflector retraction. 

6. PS-13/3/95 

For this run the nozzle-to-tip distance was set to 6" (16 d,) and 0.4” (0.5 
laterally. The ejector flow pressure was reduced to 35 psia. An essentially neutral bum 
for 1.5 seconds was achieved with a peak pressure of 360 psia. The PDPA data appeared 
to be correct. However, it had to be copied by hand and then reconstructed due to a 
computer hard drive problem. Only four valid particles were measured with a d^^ = 2.3 
|xm. Transmittance data still suggested that the ejector flow rate was too high. 

7. PS-13/4/95 

This run had a nozzle to tip distance of 6" (16 d^) and was offset 0.6” (1.6 d^ 
laterally. The ejector flow pressure was reduced to 25 psia. Due to rainy weather the 
motor took a while to light-off, but did achieve a neutral bum for two seconds near 375 
psia. The PDPA measured only one valid particle at 6.4 ^im, and the transmittance data 
looked good- 

8. PS-13/5/95 

In order to conduct a background radiation check and try an open-air, three- 
wavelength measurement, the probe with side and vertical tubes as well as the PDPA 
were not used. The motor was set up for a 6" plume length and a 0.75” offset from the 
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light beam emanating from the lab table. The neutral density filters and diode box 
pinhole had to be adjusted to prevent the diodes from saturating with the increased 
radiance. During the video replay a faint glow of the beam was visible where it came out 
of the table prior to ignition. Once the motor fired the beam was very visible as it was 
scattered by the particles in the plume. The plume appeared to be approximately 3" 
across at the location of the beam. The surprising result was plume emission at the UV 
wavelength \ = 0.313 urn, and potentially at the other wavelengths. 

9. PS-13/6/9S 

Due to the results of the previous run a light chopper was added to the setup 
between the light source and mirror and run at a speed to produce an approximately 
40 Hz square wave of alternating light and darkness. This frequency was selected 
because with the sampling rate of 300 Hz it was the highest frequency that would provide 
two to three sample points near the top and bottom, corresponding to light and no light, of 
the wave. The motor was positioned 6" (16t/J from the probe tip on centerline. The 
bum was fairly regressive with a peak of 500 psia. The PDPA had no valid particles of 
192 attempts, however the tip melted at some point during the run. Transmittance data 
indicated that, inside the probe plume, emission was not a problem since throughout the 
run accurate "no light” samples were recorded at the bottom of the square wave. 

10. PS-13/7/95 

The final firing used only the Malvern particle sizer, The Malvern was set up in a 
separate test cell with a 100 mm Fourier transform lens, providing a vignetting distance of 
approximately 5" in front of the receiver lens. A laser bandpass filter was used to 
eliminate background radiation. The motor was mounted for a nozzle-to-beam distance 
of 6" with the plume centerline 2.5" from the receiver lens. A separate LABTECH 
Notebook program was used to trigger the Malvern 1 second after 100 psia chamber 
pressure. The Malvern was initially aligned and adjusted for as low a background as 
possible. Then it was set in a program mode, awaiting external trigger to take 50 sweeps 
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in 0.4 seconds, followed by a 10 second delay, then another 50 sweep sample in 0.4 
seconds. This latter sample was to check the condition of the windows under post-firing 
conditions. The video camera was also set up looking down on the plume. The motor 
did not have much of a neutral burn and reached a peak of about 475 psia. The Malvern 
measurement was successful and occurred near the peak pressure. Post-firing inspection 
of the receiver window showed that it was partially obscured by some exhaust particle 
residue. The video replay provided an excellent opportunity to measure the plume 
dimensions that were lost on the first motor firing with PS-1. The highly visible portion 
of the plume appeared to be 12" long and 1.5" wide with the widest point at the 6" 
longitudinal position. The outer smoky portion of the plume grew from about the 3" aft 
position to a maximum diameter of about 4" near the 6" aft position, and maintained this 
width to the edge of the field of view. 
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V. RESULTS & DISCUSSION 


This chapter examines the results obtained with the Maivem. PDPA, and the 
three-wavelength apparatus. The discussion refers to the PS-1 firings on 2/25/95 and 
later. Calculations are presented to correlate the data by different means. Table 4 
presents the PS-1 firings that were used for data reduction (i.e., the first firing has not 
been included). The d, = 0.239" nozzle which was used resulted in an average chamber 
pressure (P ^) between 290 and 480 psia. This yielded underexpanded nozzle flow (exit 
pressure between 30 and 40 psia as verified with the MICROPEP equilibrium, adiabatic 
combustion/performance code run in the off-design mode for the appropriate expansion 
ratio). 


Firing 

Positiont 


3-X Ext. 

Maivem 

SEM 

Ave. 

(psia) 

Comments 

m3 



Yes 

No 


380 





No 

No 


290 


1 3/1/95 

4'70’' 

No 

Yes 

No 


375 

Melted tip 

1 3/2/95 



Yes 

No 


360 

Melted tip | 

3/3/95 



Yes 

No 


325 


3/4/95 



Yes 

No 


370 


3/5/95 

6''/0.75" 

No 

Yes 

No 

No 

400 

Bare light 

3/6/95 

6'70" 

Yes* 

Yes 

No 

No 

480 

Light chopper, 
melted tip 

3/7/95 

6''/NA 

No 

No 

Yes 

No 

470 



Table 4. Summary of PS-1 Firings Used for Data Reduction. (* - An attempt was made 
but there were no valid particles, t - Axial location/radial displacement) 


A. MOTOR 

As discussed in the previous chapter the propellant always burned in a smooth 
fashion. Periodic puffing of the exhaust from shedding of AljO, slag was not observed. 
Formerly used propellants have shown a propensity for this shedding phenomenon and 
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with such a high aluminum content in PS-1 it was anticipated. Chamber pressures varied 
between 250 and 500 psia compared to the near 300 psia calculated for ideal conditions, 
due primarily to a reduction in throat area from Al^O, deposits. 

B. EXPECTED RESULTS 

Previous experimentation had demonstrated the individual usefulness and 
limitations of the PDPA, the Malvern, and the multiple-wavelength extinction techniques 
for the solid rocket motor environment. The Malvern takes an across-the-plume 
measurement of "large” particles with d > 0.5 nm and provides number and volume 
distributions assuming there are no particles smaller than 0.5 pm. The PDPA counts 
individual particles with d > 0.5 pm passing through its probe volume and provides 
number and velocity distributions as well as representative mean diameters. It has been 
demonstrated to accurately measure the distribution of particles with d > 0.5 p,m 
suspended in liquid. Multiple-wavelength extinction measurements have also been used 
before at the laboratory. A Mie code and an assumed distribution shape (e.g., log-normal) 
permitted the particle index of refraction and distribution parameters for "small'' particles 
in the plume edges to be found. Since past work had found the larger particles only 
toward the center of the plume and the smaller particles throughout the plume, the goal of 
this investigation was to use the PDPA to obtain "large" particle distributions at various 
radial positions across the plume and to use the extinction technique at the plume edge to 
get distribution information for the ''small" particles which then could be applied across 
the whole plume. The Malvern was to be used to provide a check on the combined radial 
PDPA measurements. In order to accomplish this objective experiments were conducted 
using the probe at radial positions of 0", 0.2", 0.4", 0.6", 0.75", and 1.2" from the plume 
centerline and a final measurement was made across the plume. 

To investigate the potential for taking the obtained log-normal distribution, found 
for the small particles using the extinction technique, and curve-fitting it to a PDPA 
distribution, a spreadsheet was created in LOTUS* 1-2-3* Release 4 for Windows™. 



The PDPA distribution was manually entered for all the possible bin diameters and bins 
for additional diameters smaller than 0.5 pm were inserted for use with the log-normal 
distribution. The equation for the cumulative log-normal distribution. Equation (2.26), 
was used to determine the number of particles in each individual bin to create a number 
distribution in histogram form for comparison to that of the PDPA. An overlapping 
portion of the two distributions for diameters greater than 0.5 urn was used for the 
curve-fit. A Sum of the Squared Deviations (SSD) approach was used to obtain a single 
number to evaluate the curve-fit. The SOLVER routine in 1-2-3 was used to adjust d „ 
Gg, and N to minimize the SSD. Unfortunately, even with constraints the routine could 
not accomplish this. It could, however, work with a single parameter if the other two 
were fixed. By manually adjusting the parameters the problem was readily apparent. 
Many log-normal distributions can provide an equally good fit. Thus, curve-fitting will 
not result in a unique solution. 

The Mie code FORTRAN program at NPS, obtained from Cashdollar [Ref. 12], 
was formerly only available on the main-frame computer. Therefore, to speed up the data 
reduction and analysis, it was copied and modified for use on a personal computer with 
Microsoft MS-Fortran 4.1. To simulate experimental conditions in the 1-2-3 spreadsheet 
and examine expected transmittance values for different distributions and at different 
wavelengths "miescat3.exe’' (see Appendix C) was used and a macro was written for 
1-2-3 to be able to properly import the output file (see Appendix D). What was more 
important was also to be able to examine the InT or Q ratios as the mass of "small" 
panicles was varied with respect to the mass of the "large" panicles. Finally, for data 
reduction (i.e., finding the index of refraction m, and that best fits the 
experimentally derived transmittances) the program "miescatb.exe" was developed (see 
Appendix E). 

The most optimistic result desired was to be able to take the end-product, 
experimentally-derived distributions for various radial positions across the plume and 
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conduct a mass calculation for comparison to the calculated mass that should be present 
from burning the propellant. 

C. MALVERN 

The Malvern data were typical of past observations with propellants of various 
composition. Beam steering occurred as evidenced by the increased intensities on the 
first six detector rings. The voltages on these rings were deleted before the size 
distribution was calculated. Elimination if these inner six rings does not affect the 
distribution when all the particles have diameters less than the expected value of 
approximately 30 |um. In agreement with previous data obtained with the subscale motor 
[Ref. 18], the Malvern measured -25% of the volume in particles with d < 2.0 [im, the 
remaining -75% of the volume in particles with if = 4 - 5 (xm, and a ~ 6 nm. An 
obscuration of 86% (therefore 7 = 0.14) was recorded, with beam-steering present, and a 
(fjj = 2.59 |xm was determined. 

In order to conduct a check on the above value of T, Dobbins' formula (Equation 
2.30) was used. First, a stagnation enthalpy energy balance equation was used to estimate 
the ratio of the mass flow rates of air to propellant gas in the plume at the measurement 
location for use in calculating the mass concentration of particles. 


+ »iairCp_,r,^, — 


(5.1) 


where 


Cp^, 

T,., 

Cp^ 


= mass flow rate of propellant gas 

= specific heat for propellant gas = 1431 J/kg-K from MICROPEP 
= stagnation temperature of the propellant gas = 3584 K from 
MICROPEP 

= mass flow rate of entrained air 

= specific heat for entrained atmospheric air = 1004 J/kg-K 
= stagnation temperature of entrained air = 520 K 
= specific heat of the mixture of gases = 1217 J/kg-K from 
average of gas and air 

= stagnation temperature of the mixture of gases-I500K [Ref. 19 
at 16<fj 
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Dividing by the mass flow rate of propellant gas, the ratio of the mass flow rates was 
calculated- The mass flow rate of entrained air was approximately 2.5 times that of the 
propellant gas. Therefore, the mass flow rate of the naixture of gases was 3.5 times the 
mass flow rate of the propellant gas alone. The equation used for C„ was 


Cm 



(5.2) 


where 

Vmix'^ volume flow rate of the mixture of gases 
density of the mixture of gases 
Pm. - pressure of the mixture of gases - 14.7 psia 
R„^= gas constant for the mixture of gases ~ 54 ft-lbf/lbm-°R 
temperature of the mixture of gases - 1500K 

First, the mass flow rate of the propellant had to be calculated using 

mp = ppA*ap" (5.3) 

where 

rkp = mass flow rate of the propellant 

= propellant density = 0.06665 Ibm/in^ from MICROPEP 
Ap s= propellant burning surface area = 4.8 in* average during steady bum 
a = propellant burning rate pressure coefficient 
= chamber pressure = 400 psia 
n = propellant burning rate pressure exponent 

Next, the mass flow rate of AljO, was found by using the mass percentage of aluminum 
in the propellant (20%) times the ratio (1.89) of the molecular weights of Al^O, (101.96 
g/mol) to aluminum (2 x 26.98 g/mol) times the mass flow rate of the propellant found 
using Equation (5.3). Finally all values were plugged into Equation (5.2) to find C„= 1.4 
X 10* Ibm/in’. Using Equation (2.30) with L = 3", fi = 2. p = Pai 2 oj= 0.138 Ibm/in* at 
1500 K, and = 2.59 pm a transmittance of r= 0.40 was calculated. This is much 
higher than the Malvern measured value of 0.14, most likely due to beam steering as 
discussed above. Beam steering is the refraction of the beam due to the thermal gradient 
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in the path of the beam. This causes the central detector, which measures the obscuration, 
to receive less intensity than it should and the innermost detector rings to receive too 
large an intensity. 

Another calculation was performed to compare the mass in the observed modes 
for the calculated 4,^ of 2.59 pm. For the Malvern volume distribution two modes were 
apparent, near 1 and 4.5 pm. Therefore, it was assumed that all the particles were either 1 
or 4.5 pm and the number ratio of 1 to 4,5 pm particles was found to be 24.3 using: 

N, d32(4.5)^-(4.5)^ 

N 4.5 (l)5-rf32(l)' 


where N, is the number of 1 pm particles and N,, is the number of 4,5 pm particles. 
Then, the mass percent of the 1 pm particles to the total mass of the particles was 
calculated as 21.1 % using 


Vi 

Vtouil 


(5.5) 


where V, is the volume of 1 pm particles and is the volume of all the particles. Thus, 
the Malvern data implied that most of the mass was contained in particles with diameters 
greater than I pm. This is evident in Figure 10 of Appendix G. 

D. PDPA 

Six attempts were made to obtain measurements at various radial positions in the 
plume with the PDPA, from the centerline out to 1.2" (3 d,). Both centerline firings 
resulted in no valid particles. On both of these runs as well as the 0.2" offset run the 
probe tip was melted. Therefore, flow deflection and/or a strong detached normal shock 
could have caused the lack of "large" particles in the probe by breaking them up into 
"small" particles. 
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The Malvern number and volume distributions were compared to an equivalent 
across-the-piume distribution formed by scaling the 0.2", 0.4", and 0.6" offset PDPA 
distributions- This scaling was done by assuming the plume was axisymmetric and could 
be divided up into concentric zones. Within each zone the distribution was assumed to be 
constant. Since there were three PDPA distributions that could be used three zones were 
formed with the constraint that each zone could contain only one PDPA distribution. 

This forced the inner zone to be constrained to an outside radius of 0.2 - 0.4". Similarly, 
the middle zone extended from the outside of the inner zone to 0.4 - 0.6" and the outer 
zone extended from the outside of the middle zone to 0.6 - 0.75". Each distribution was 
then multiplied by twice the difference between its zone's inside and outside radii. This 
was the scale factor. Next, the scaled distributions were combined and the total number 
of particles found. This total was then divided into each number value and multiplied by 
100 to obtain number % at a given diameter. Finally, the number %’s were added 
together according to the bin diameters used by the Malvern and plotted with the Malvern 
distribution. An additional step was required for the volume distribution. Each scaled 
number was multiplied by (P. Then when combined, the total volume of particles was 
found rather than number. The PDPA and Malvern plots of the number and volume 
distributions (see Figures 10 -13 in Appendix G) agreed that there was a mode at < 1 
and another mode between 4 and 5 pm. The PDPA volume distribution highlights 
individual larger particles that it measured, whereas individual particles are essentially 
lost by the Malvern as it looks at an ensemble of particles. 

The PDPA distributions also showed that there were almost no "large" particles (el 
> 0.5 pm) outside of 0.6" off the centerline. This was in agreement with the visual 
appearance of the plume. The brighter plume core area appeared to have approximately a 
1" diameter. 

There are several ways to improve the future performance of the PDPA when used 
in the plume environment- First and foremost, if the probe tip were cooled or made of a 
more thermally resistant material the tip erosion problem could be minimized. A shorter 
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time for tip exposure would also help in this regard. The ejector flow needs to be 
optimized to prevent sucking in too many of the small panicles since this biases the data 
toward smaller panicles. It would be very difficult to obtain isokinetic sampling because 
of the required shon explosive time. The probe needs to be redesigned with a smaller 
frontal area, probably circular, and with a diameter no larger than the present maximum 
tip diameter. This would cut down on the plume obstruction and remove some of the 
excess internal volume that presently causes recirculation within the probe. If fiber optics 
were used, the large side tubes could also be removed, funher improving the plume 
obstruction problem. With subscale motors it may also be possible to not use the probe, 
allowing the PDPA to traverse the plume cross-section. 

E. THREE-WAVELENGTH EXTINCTION 

Transmittances were calculated for seven firings and used to obtain "best fit" 
values of index of refraction m, and for a log-normal distribution. The following 
formula was used: 


T- 


Vp-Vo 
V/-Vo 


(5.6) 


where V is the diode output in volts for particle obscured light, Vg is the diode output for 
no light, and V, is the diode output for unobstructed full light. The no-light reading was 
made prior to firing by covering the hole in the bottom of the lab table. The full light 
value was found by taking an average of the first two seconds of the run (i.e., before 
ejector gas activation which typically reduced the diode reading slightly). The particle 
obstructed reading was taken over a flat minimum during steady-state burning. 

Once the three transmittances were calculated for a firing they were inserted into 
the data input file ’'miescat6.dat” and iteration intervals were chosen for m and a . The 
data reduction file, "miescat6.exe", was run and provided the best t/jj for each 
combination of m and as well as a calculation of the degrec-of-fit. This value was the 
Sum of the Squared Deviations (SSD) between InT ratios and Q ratios. The closer SSD 
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approached zero the better the fit. Table 5 presents all of this information. The SSD 
values are presented as integers rather than the eight digits to the right of the decimal that 
they really were (i.e., 237 is really 237E-08). The 3/5/95 firing was different from the 
others because the probe was not used and the plume was definitely emitting at 0.313 nm. 
Therefore, a transmittance value was not calculated at this wavelength and a fit could not 
be obtained. 

It was discovered that the process used to iterate toward the "best-fit” values was 
convergent for o and dy^. but not for m. Therefore, since a reasonable value for m of 
AJjOj in the plume is 1.70 - Oi, this value was used as a starting point to find a reasonable 


Firing 

Position 

T-in 

7,44 




dy, 

SSD 

2/25/95 

6'7 1.2" 

0.0557 

0.0533 

0.0958 

1.70-0i 

1.82 

0.6 

237 

3/1/95 

4'’/0" 

0.2891 

0.3304 

0.4905 

1.70-0i 

1.7 

0,44 

1,708 

3/2/95 

6"/0.2" 

0.4957 

0.5054 

0.5666 

1.70-0i 

2.1 

0.69 

925 

3/3/95 

6"/0.4" 

0.3026 

0.425 

0.5567 

1.70-0i 

3.24 

0,36 

292 

3/4/95 

6'7 0.6’' 

0.2518 

0.3315 

0.4316 

1.70-01 

3.38 

0.58 

253 

3/5/95 

6'7 0.75" 

N/A 

0.794 

0.7184 

N/A 

N/A 

N/A 

N/A 

3/6/95 

6'7 0'’ 

0.1711 

0,2268 

0.3007 

1.70-0i 

3.9 

0.84 

207 


Table 5. Three-Wavelength Extinction Results. 


value for and djj. Then, if a better fit was desired, m had to be varied over its entire 
possible range in 0.01 increments as was stepped in its smallest increment of 0.01 
about the reasonable value and the minimum SSD set chosen. This entire somewhat 
lengthy process was not done, however some excellent fits were found. For the 2/25/95 
firing m = 1.79 - 01,0^ = 1.84, dj2 = 0.55 resulted in an SSD= lE-08. For the 3/4/95 
firing m= 1.85 - Oi, = 3.38, djj = 0.50 resulted in an SSD = 2E-08. The large values of 
and djj relative to the 0.5 |jm particles is really indicative of a significant amount of 
the mass in "large" particles. It should be noted that only two independent ratios were 
available using three wavelengths. Kim et al. [Ref. 20] used six wavelengths for five 
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independent ratios and determined a ~ 1.50 and rfjj- 0.150 pm in the edge of the 
plume. Thus, if further work is to be done, it should use as many wavelengths as possible 
with a large difference between X and X 

An uncertainty analysis was conducted on the transmittances, then applied to a 
single transmittance and propagated through the fit process to determine how much the 
distribution parameters changed. Using Holman's [Ref. 21] approach the transmittance 
uncertainty was found with 



where Wj, w^, and are the uncertainties in the transmittance, particle obscured 
diode output, no-light diode output, and full-light diode output, respectively. Using the 
values from the 3/4/95 data set and Wy^were set equal to 0.0397, the 1 o 

deviation on the 546 diode between 0 and 2 seconds elapsed time. This resulted in a 
transmittance uncertainty of 0.008974. This was added to to create a perturbed set, 
then the fit process was conducted using m = 1.70 - Oi. This produced a = 5.85, 
tijj = 0.71, and an SSD = 8E-08. If this is compared to the original fit in Table 5 it is 
obvious that these values correspond to two unique distributions. Therefore, the multiple 
wavelength extinction technique is highly sensitive to any experimental uncertainties. 

The 3/5/95 bare-light data (no probe) were used to compare the mass of small 
particles to large particles. Since the light beam was 0.75" from the centerline of the 
plume and appeared in the video replay to be in the smoky portion of the exhaust, it was 
assumed to be in the small-particle region. Equation (2.30) was used to detemine C„. 
The plume was estimated to be 3" in diameter, so L was calculated to be 2.60". The 
distribution parameters found using the fit routine with the 2/25/95 data taken in the 
plume edge were used with the Mie code to find Q. However, the transmittance values 
were not used because the high probe purge rate resulted in a high concentration of small 
particles in the probe (note low values of T). Using the 3/5/95 rjso value, a C„ of 1.2E-07 
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Ibm/in’ was calculated. When compared to the total C„ of 1.4E-06, calculated above, 
this showed that the small particles accounted for only 8% of the total mass. Similarly, 
using the 3/5/95 value a C„ of 6.4E-08 Ibm/in^ was calculated which corresponds to a 
small particle mass of only 4% of the total. When the mass distribution is plotted for the 
small particles (see Figure 18 in Appendix G) notice that only 30% of the "small panicle" 
mass is in particles with diameter less than 0.5 pm. Therefore, the mass % of particles 
with d < 0.5 pm is less than 3%. Thus, accurate measurement of the size distribution in 
the plume should be possible using only the PDF A. This greatly simplifies the 
diagnostics required to make the needed plume panicle size measurements. 

F. SEM ANALYSIS 

Several samples of the particles coating the inside rear of the probe were taken for 
analysis under a scanning electron microscope (SEM). First, the panicles were 
sedimented in a beaker with acetone to try to remove any other matter. Next, they were 
scraped from the bottom and sprinkled on carbon painted mounts. Finally, the mounts 
were placed in a vacuum container where carbon fibers were vaporized to carbon coat the 
samples. This was to help with the sample charging problem. The samples were 
examined qualitatively under the SEM for an idea of the general size distribution. See the 
photographs in Figure 19 of Appendix G. It was noted that there were many 1 to 2 and 4 
to 6 pm particles and few smaller or larger particles. This is consistent with the measured 
distribution, using the Malvern and the PDPA, as well as the various mass calculations. 
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VI. CONCLUSIONS 


The PDPA and MaJvern distributions agreed in the observed modes near 1 and 
4.5 nm. SEM pictures of collected particles were in good qualitative agreement with the 
Malvern Sauter mean diameter = 2.59 pm. Data analysis indicates that less than 3% 
of the total mass of the particles was contained in particles with diameter d < 0.5 pm. 

This means that the PDPA, which can typically measure particles down to a minimum 
diameter of 0.5 pm with a dynamic range id^'.d^^ of 50:1, can be used by itself to 
determine the particle size distribution. Although multiple-wavelength measurements 
were not required, they were found to be very sensitive to inaccuracies in the measured 
transmittances. The modified particle collection probe allowed simultaneous PDPA and 
multiple-wavelength measurements. However, the design should be improved in several 
ways. The frontal area of the probe should be reduced to lessen plume obstruction. One 
option is to use a small cylindrical cross-section and use fiber optics with the PDPA, 
thereby removing the large side tubes which also obstruct the plume. In addition, the tip 
exposure time should be reduced or a more thermally resistant material chosen for the Up. 
Further studies should be conducted to check shock formation in and around the tip. A 
"rake" of PDPA miniature fiber optic probes may be a reasonable design to investigate for 
taking simultaneous radial measurements. Traversing such a rake across the plume is 
another option. 
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APPENDIX A. MICROPEP OUTPUT 


The following pages contain a printout from the equilibrium, adiabatic 
combustion/performance code MICROPEP. This particular run was for propellant PS-1 
with a chamber pressure of 400 psia and a nozzle area ratio of 2.33 operating off-design 
for an ambient sea level pressure of 14,7 psia. The area ratio pertains to the d, = 0.239" 
nozzle used for most of the experiments with PS-1. Note that the exit pressure is about 
40 psia, indicating underexpanded flow. 
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Feb. 1994 - Modified by A. HcAtee 

Naval Postgraduate School, Monterey, 
NEWPEP »• 

• PS-1 Ae/At=2.33 » 

INGREDIENT MASS 

(gm) 

AMMONIUM PERCHLORATE (AP) 70.00 

ALUMINUM (PURE CRYSTALINE) 20.00 

HTPB (SINCLAIR) 10.00 


COMPOSITION 


DENSITY 
(kg/li.**3) 

1948.6650 ICL 4H 

2701.5590 lAL 

918.9728 103H 73C 


0 13.0 

VOLUME PERCENT OF INGREDIENTS (IN ORDER) 

66.2684 13.6572 20.0744 

THE PROPELLANT DENSITY IS .06665 LB/CU-IN OR 1.8448 GM/CC 
THE EQUIVALENCE RATIO IS 1.3260 

NUMBER OF GRAM ATOMS OP EACH ELEMENT PRESENT IN INGREDIENTS 




2.393073 C 


• >*«***********>*>***««*«***CfU^Q£I{ RESULTS FOLLOW’ 

TEMP 

(K) 

3584.' 

DAMPED AND UNDAMPED SPEED OF SOUND= 848.508 AND 1063.794 n/sec 

54036.950 J/kmol-l 


PRESSURE 
(HPa/ATM/PSI) 
758/ 27.22/ 400 


ENTHALPY ENTROPY CP/CV SGAMHA Pi/ni 

()tJ/kg) (kJ/)tg-K) (MPa/kmol) 

•1757.6980 9.398 1.1620 1.1212 84365.270 


“liquid,S=solid) 
.84621 H2 
.33919 A1203* 
.06428 HO 
90E-03 0 
22E-03 AlO 
45E-04 A120 
25E-05 OCl 
18E-05 A102 
81E-05 NH3 
17E-06 NHO 


.68459 CO 
.29499 N2 
.04771 C02 
82E-03 A1C12 
99E-03 AlHO 
80E-04 A1C13 
47E-05 A1202 
07E-05 HOCl 
81E-05 H02 
92E-06 AlHO 


.48272 H20 
.22731 H 
.03125 AlCl 


93E-05 COCl 
75E-05 HH2 
68E-06 CH20 

27.721 gm/mole 


.45869 HCl 
.08121 Cl 
.01142 AlOCl 
92E-03 A1H02 
46E-03 A1 


12E-a 


C12 


99E-05 CHO 
82E-05 NH 
93E-06 CNH 


PRESSURE 
(MPa/ATM/PSI) 
.276/ 2.72/ 40 


•‘’►EXHAUST RESULTS FOLLOW*******“““““* 

ENTHALPY ENTROPY CP/CV SGAMHA 

(kJ/kg) (kJ/kg-K) 

.02 -3746.0840 9.398 1.1556 1.1201 


DAMPED AND UNDAMPED SPEED OF S0UND= 732.222 AND 938.454 n/sec 
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SPECIFIC HEAT (MOLAfi) OF GAS AND TOTAL=39437.610 55329.410 J/kmol-K 

NUMBER MOLS GAS AND CONDENSED= 3.1425 .3642 


(*=liquid,i=solid) 


.8B420 H2 
.36417 A1203* 
5.05E-02 Cl 
1.82E-03 0 
4.60E-04 AlHO 
9.42E-05 A1 
6.01E-06 CHO 
3.41E-06 COCl 
1.15E-06 A102 


.67999 CO 
.29725 N2 
2.15E-02 HO 
1.27E-03 A1C12 
4.01E-04 02 
3.24E-05 C12 
5.81E-06 AlH 
2.96E-06 NH3 
1.09E-06 H02 


.53227 HCl 
.12299 H 
..B5E-03 AlCl 
..20E-03 NO 
:.01E-04 AlO 
..S3E-05 A120 
;.79E-06 HOCl 
..86E-06 A1202 
..03E-06 NH 


.48500 H20 
.05240 C02 
3.03E-03 AlOCl 
8.34E-04 A1H02 
1.33E-04 A1C13 
B.73E-06 OCl 
4.70E-06 N 
1.42E-06 NH2 
l.OOE-06 CNO 


THE MOLECULAR WEIGHT OF THE MIXTURE IS 28.517 gm/mole 
THE GAS CONSTANT IS 291.56 J/kg-K 

TOTAL HEAT CONTENT (298 REF) = 4788.339 kJ/kg 

SENSIBLE HEAT CONTENT (298 REF)“ 4456.853 kJ/kg 


An exact nettiod for determining throat conditions was used 

The frozen & shifting STATE gammas for the throat are: 1.1600 1.1199 

GAMMA NU shown below is the gamma for the chamber to throat PROCESS. 


.»**»**.**PERF0RMANCE: FROZEN ON FIRST LINE, SHIFTING ON SECOND LINE.******* 
.***.***»«.»*.SHIFTING TO THROAT/FFOZEN AFTERWARDS ON THIRD LINE************* 


SPECIFIC GAMMA T* P* C* ISP* 
IMPULSE KU 

(sec) (K) [MPa) (m/s) (sec) 

222.8 1.1630 3314. 1.577 1539.2 

226.8 1.1202 3416. 1.601 1560.4 195.B 

224.9 1.1202 3416. 1.601 1560.4 195.8 


Pe D-ISP A*/ni Te Cf 
(gm-s/ (cn**2/ 

(MPa) cin**3) kg/s) (K) 

.276 10192.710 5.5805 2590. 1.2718 

.297 10238.540 5.6571 2937. 1.2602 

.297 10140.370 5.6571 2702. 1.2482 


FROZEN & SHIFTING KINETIC ENERGY OF EXHAUST 674708. 781129. 
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APPENDIX C. MIECODE"MIESCAT3.FOR" 


On the following pages is the Fortran program to calculate extinction coefficients 
Q, average extinction coefficients Q, and ratios of average extinction coefficients for the 
inputs given in "miescatS.dat". The original program was run on the main-frame 
computer at NPS. Since the Combustion Research Lab (CRL) did not have a network 
connection to the main-frame the program was copied and modified to mn on an IBM 
compatible PC using the Microsoft* MS® Fortran 4.1 compiler at the CRL. 

The data block in the original program followed the end subroutine. Since the 
source code had to be compiled and linked prior to execution, it was inconvenient to have 
to repeat these steps each time the inputs were changed. Therefore, the data block was 
removed and put in the separate file "miescat3.dat'' and the read statements modified in 
the source code. 
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BY KENNETH CASHDOLLAH, 1976 
REVISED BY THOMAS WELDON. JULY, 1977 
BASED ON PROGRAM MICOER BY 0 D LITTON AND Z J FINK 
REVISED BY LT JOHN MANSER, USN 1994-1995 (UPDATED FOR 
MS-FORTRAN 4.1,DATA BLOCK WITH WAVES PUT INTO 
EXTERNAL FILE) 


"MIESCAT3.EXE" CALCULATES EXTINCTION COEFFICIENTS AND RATIOS FOR 
LIGHT TRANSMISSION THROUGH DUST PARTICLES WITH COMPLEX REFRACTIVE 
INDEX AND A LOG-NORMAL SIZE DISTRIBUTION 


"MIESCATl-DAT" INPUT DATA FILE FORMAT 


VARIABLE NAMES AS IN PROGRAM 

ROW 1.WAV(l) ,WAVt2) ,WAV(3) 

ROW 2.NDATA 

ROW 3.RDUST.RMED 

ROW 4.L2,L3,DS.DELDS 

ROW 5.ND,SIGMA,SD32,SDL32 

FORMATS AND VARIABLE DESCRIPTION 

ROW 1_F6.4 (WAVELENGTH 1),2X,F6.4 (WAVELENGTH 2).2X, 

F6.4 (WAVELENGTH 3) 

ROW 2_13 (DATA SET IDENTIFIER, IF EQUAL TO 999 SUBSEQUENT 

DATA IS NOT PROCESSED). 

ROW 3....2F6.3 (COMPLEX REFRACTIVE INDEX OF DUST PARTICLES), 
2X,F6.3 (REFRACTIVE INDEX OF THE MEDIUM). 

ROW 4....13 (L2, THE NUMBER OF INCREMENTS TAKEN AT A SET DELTA 
)(THE NUMBER OF TIMES DELTA IS DOUBLED. L3),F6.3 
(STAR DIAMETER),F6.3tINITIAL VALUE OF DELTA). 

ROW 5_I3(NUMBER OF INCREMENTS TAKEN FOR LOG-NORM), F4.1 (SIGM 

F5.2(STARTING DIAMETER).F5.2(DELTA). 


COMPLEX B,RDUST,DI,DELDI.RX,RINDX,DIS,DELIS 
COMPLEX COEFR.DPRdOOO) , XJR( 10001 ,DJR(lOOO) ,PR(:000) 
DIMENSION XY(IOOO),WAV(3) 

DIMENSION DP(1000),XK(1000).XJ(1000),DJ(1000),P(1000) 
DIMENSION RATl(lOO).RAT2(100).RAT3(100) 

DIMENSION PLOTDY(IOO) 

COMPLEX A(IOOO),B(1000),0(1000) 

COMPLEX ALPH,XH(1000),DH(1000),E(1000),DE(1000) 
COMMON/ADQ/DQ(500,2) 

COMMON ALLQ(100,4,3) 

COMMON/KPLT/RINDX(3),WAVLt3) 


66 









COMMON DIAM32(100) 


IPEN = 4 

C CALL SHERPAl'FILENAME','M',IPEN) 

OPEN (UNIT=1, FILE='HrESCAT3.DAT') 
READ (1,103)WAV(1),WAV(2),WAV(3) 
103 FORMAT!F6.4,2X,F6.4,2X,FS.4) 

DO 8 NDATA=I,10 
READ (1,105)NDAT 
105 FORMAT!13) 

IF(NDAT.EQ.999) GO TO 9 


C RDUST IS COMPLEX REFRACTIVE INDEX OF DUST PARTICLES, WITH MINUS 
C RHED IS REAL REFRACTIVE INDEX OF MEDIUM 
READ (1,32)RDUST,RMED 
32 FORMAT(2F6.3,2X,F6.3) 

C RX IS RELATIVE REFRACTIVE INDEX, (N-IK)/H 

RX=RDUST/RMED 

C IMAGINARV PART OF R CHANGED TO POSITIVE TO CONFORM TO EXPANSION 
C FORMULA IN PROGRAM 

R=CONJG(RX) 

R1NDX(1)=RX 

RINDX(2)=RX 

RINDX(3)=RX 

C L2*L3 = NUMBER OF SIZE INCREMENTS 

C INITIAL DIAMETER IS DS, INITIAL VALUE OF DELTA IS DELDS 
C 

READ (1,34)L2,L3,DS,DELDS 
34 FORMAT{2I3,2F6.3) 

C 

C INITIAL VALUE OF COMPLEX DIAMETER IS DIS 
DIS=CHPLX(0.0,DS) 

DELIS=CHPLX(0.0,DELDS) 



READ (1,I2)ND,SIGHA,SD32,SDL32 
12 FORMAT(I3,F5.2,2F5.2) 

HAVE IS WAVELENGTH OF LIGHT IN 
WAV(l)-0.3650 
WAV(2)=0.5461 
WAV(3)=0.8520 
INITIAL VALUES FOR ARRAY 
DO 51 K=l,3 
DO 52 J=l,2 
ALLQd , J,K)=0.0 
ALLQ(1,J+2,K)=0.0 
DO 53 1= 2,100 
ALLQd,J,K)=100.0 
ALLQCI,J+2,K)=l.0 
53 CONTINUE 
52 CONTINUE 
51 CONTINUE 

DO 54 1=1,100 
DIAM32(I)=0.0 
RAT1(I)=0.0 


MICRONS 
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FUiT2(I)=0.0 
RAT3(I)=0.0 
PLOTDY(Il-0.0 
54 CONTINUE 
DO 4 HI=1,3 

PRINT 106,NDATA 

106 PORMAT(/' DATA SET NUMBER',13) 

PRINT 35 

35 70RMAT(' HIE COEFFICIENTS FOR HONODISPERSE DUST PARTICLES') 

PRINT 41,RDUST,RHED,RX 

41 FORMAT(' RDUST='.2F7.4,/' RMED=',F7.4,/' REL. REF. INDEX-',2F7.4) 
PRINT 104,R 

104 FORMATC' R=',2F7.4) 

WAVE=WAV(NI) 

PRINT 30,WAVE 

30 FORMAT!' WAVELENGTH F7.4,2X,'MICRONS'/) 

PRINT 25 

25 FORMAT(SX,'A',12X.'D',13X,'QABS',IIX,'QEXT'/) 
WAVL(NI)=WAVE*10000.0 


C INITIAL VALUES SET 

C 

C MAX VALUE FOR D IS ( 1000*WAVE)/(1.3»PI*RMED) = 244*WAVE/RMED 
C D IS REAL DIAMETER IN MICRONS,DELD IS INCREMENT 
C DI=1*D IS COMPLEX DIAMETER, DELDI IS INCREMENT 
C X IS REAL SIZE PARAMETER X= PI»D1AMETER*RMED/WAVELENGTH 
C ALPH IS COMPLEX SIZE PARAMETER, ALPH-I*X 

D=DS 

DELD-DELDS 

DI=DIS 

DELDI-DELIS 


M-1 

DO 2 1=1|l2 

X-3.14159*D/MAVE*RMED 

ALPH»3.14159*DI/WAVE*RHED 

XJ(1)= SIN(X)/X 

XJ(2)- SIN(X)/(X»*2)- COS(X)/X 

XJR(1)-CSIN(R*X)/(R*X) 

XJR(2)=CSIH(R*X)/((R»X)*»2)-CC0S(R»X)/(R*X) 

XH (1)=CEXP(ALPH)/ALPH 

XH(2)-((-CEXP(ALPH))/X)-((0..1.)*CEXP(ALPH)/(X**2) ) 
DJ(1)= COS(X)/X- SIN(X)/(X»*2) 
DJR(l)=CCOS(R*X)/(R*X)-CSIN(R*X)/(CR*X)**2) 
DH(1)=CEXP(ALPH)/X+((0.,1.)»CEXP(ALPH)/(X**2)) 
QEXT=0.0 
QSCA-0.0 
C 

C LI = ITERATIONS FOR SPHERICAL FUNCTIONS. L1=1.3*X 
L1=1.3»X 

IF(Ll.LE.lO) Ll=10 
DO 1 L=1,L1 
CO=(FLOAT(L))/(2»L+l) 

COE-(FLOAT(L+1] )/(2»L-M) 
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C0EF=1FL0AT(2*L+I))/X 
C0EFR=(FL0AT(2*L+1))/(R*X) 

XJ(L+2)=COEF*XJ(L+1)-XJ(L) 

XJR(L+2)=C0EFR*XJR(L+1)-XJR(L) 

XHCL*2)-»C0EF*XH(L+1)-XH1L) 

DJ(L+1)=C0*XJ(L)-COE*XJ(L*2) 

DJR(L+l)=C0*XJR(L)-C0E*XJR(L-^2) 

DH( L-t-1) =CO*XH (L)-COE*XH( L+2) 

P(L)-X*XJ(L+1) 

PR(L)=(B*X)*XJR(L+1) 

DP(L)=XJCL+1)+X*DJ(L+1) 

DPR(L)=XJR(L+1)+(R*X)*DJR(L+1) 

E(L1=X*XH(L+1) 

DE(L)=XH(L+1)»X*DH(L+1) 

A(L)=((DPR(L)*P(L))-(R*PR(L)»DP(L)))/((DPR(L)*E{L))-(R*PR(L)*DE(L 
1))1 

B(L)=((R*DPR(L)»P(L))-(PR(L)*DP(L)))/((R*DPR(L)»E(L))-(PR(L)*DE(L 
1))) 

C(L)=A(L)+B(L) 

GAHMA=CAaS{A(L)) 

BETA=CABS(B(L)) 

XK(L)=(FLOAT(2*L+l))*(2./(X*»2))*((GAHMA**2)+(BETA**2)) 

Xy(L)=‘( FLOAT! 2*L*1) )*(2./(X**2) )*(REAL(C(L) ) ) 

QEXT=QEXT+XV(L) 

QSCA=QSCA+XK(L) 

1 CONTINUE 
QABS^QEXT-QSCA 
DQ(H,1)=D 

DQ(H,2)=QEXT 

PRINT 22,X,DQ(M,1).QABS,DQ(M,2) 

22 FORMAT(IX,F8.3,6X,F7.3,9X,F6.3,9X,F6.3) 

C 

C D AND QEXT FOR MONODISPERSE PUT INTO PLOTTING ARRAY 
C EVERY FOURTH POINT CALCULATED IS PUT INTO ARRAY- ALLQ 

XH=M 

AREL=XK/4.0 

MINT=M/4 

XMINT=HINT 

IF (AREL.NE.XMINT) GO TO 60 
IF (D.GT.1.20) GO TO 60 

C NUMBER OF POINTS FOR PLOTS 1 AND 3 IS MPTSl 
NPTS1=MINT 
ALLQ tMINT+1,1,NI)=D 
ALLQ(HINT+1,3,NI)=QEXT 
ALLQ(HINT+1,1,NI)=D 
ALLQ(HINT+1,3.NI)=QEXT 
60 CONTINUE 

C INCREMENT SIZE PARAMETERS D=DIAMETER 

D=D+DELD 
DI=DI+DELDr 
M-M+1 

2 CONTINUE 

C PRINT 36 

36 FORMAT (IH ) 

DELD=2.0*DELD 

DBLDI=2.0»DELDI 

3 CONTINUE 
PRINT 106,NDATA 
PRINT 235 
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235 FORHATf' HIE COEFFICIENTS FOR LOG-NORMAL DUST DISTRIBUTION') 
PRINT 236,SIGMA 

236 FORMAT (' SIGMA=',F4.2) 

PRINT 237.RX 

237 FORMATC RELATIVE REFRACTIVE INDEX=',2F7.4) 

PRINT 30,WAVE 

CALL QAVG(M,NI,SIGMA,NPTS2,SD32,SDL32,ND) 

K1-ND*1 

DO 499 N=1,K1 

PLOTDY(N)=ALLQ(N,4,NI) 

499 CONTINUE 

IF(NI.EQ.2) CO TO 502 
IF(NI.EQ.3) GO TO 503 

501 CALL PLOTD(DIAM32,PLOTDY,41,.FALSE.,'LINLIN','NOLEGEND', 

+ ' ','Diameter, microns','Extinction Coefficient') 


C 502 CALL PLOTD(DIAM32,PLOTDY,41,.FALSE.,'LINLIN','NOLEGEND', 

0 + ' ','Diameter, microns','Extinction Coefficient') 

c GO TO 4 

C 503 CALL PLOTD(DIAM32,PLOTDY,41,.TRUE. ,'LINLIN','NOLEGEND', 

C + ' ','Diameter, microns','Extinction Coefficient') 

4 CONTINUE 

PRINT 106,NDATA 

7 CONTINUE 
PRINT 500 

500 F0RMAT(' QBAB RATIOS FOR 3 WAVELENGTHS') 

PRINT 238,SIGMA 
PRINT 237,RX 

PRINT 510,WAV(3),WAV(1),WAV(3),WAV(2),WAV(2),WAV(1) 

510 FORMAT(/2X,' D32 ',3(5X,'QBAB(',F6.4,')/QBAB(',F6.4,')')/) 

K=ND+1 
DO 600 N=2,K 

RAT1(N)»ALLQ(N,4,3)/ALLQ(N,4,1) 

RAT2(N)=ALLQ(N,4,3)/ALLQ(N,4,2) 

RAT3(N)-ALLQ(N,4,2)/ALLQ(N,4,1) 

PRINT 650,DIAH32(N),RAT1(N),RAr2(N),RAT3(N) 

600 CONTINUE 

650 FORMATllX,F6.3,13X,F7.4,23X,F7.4,23X,F7.4) 

PRINT 350,NDATA 

350 FORMAT (/'••**••***•*****•*•*•***•** END OF DATA SET NUMBER',13/) 
0 CALL PLOTD(DIAM32(2),RAT1(2),40,.FALSE.,'LINLIN','NOLEGEND', 

C + ' ','Diameter, microns','Extinction Coefficient') 

C CALL PLOTD(DIAH32(2) ,RAT2( 2) , 40,. FALSE.,'LINLIN' , 'NOLEGEND', 

C + ' '('Diameter, microns','Extinction Coefficient') 

C CALL PLOTD(DIAM32(2),RAT3(2),40,.TRUE. ,'LINLIN','NOLEGEND', 

C + ' '('Diameter, microns','Extinction Coefficient') 

8 CONTINUE 

9 CONTINUE 

C CALL DONEPL 


SUBROUTINE QAVG(H,NI,SIGMA,NPTS2,SD32,SDL32,ND) 

C 

C BY KENNETH CASHDOLLAR, 1976 

C 

C QAVG CALCULATES D32 AND QBAR FOR A GIVEN QEXT VS D AND 
C A GIVEN LOG-NORMAL SIZE DISTRIBUTION 

C 
C 

C INPUT DATA = ND,SIGMA,D32,DEL32 
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COMMON ALLQ(100,4,3) 

COHMON/ADQ/DQ(500.2) 

COMMON DIAM32(100] 

C D0(I,1)=D, DQ(I,2)=0EXT 

C NI=WXVELENGTH OR REFRACTIVE INDEX CHANGE 
C ND IS NUMBER OF D32,QBAR CALCULATIONS 
C SIGMA IS THE GEOMETRIC STANDARD DEVIATION 
C D32 IS THE SURFACE WEIGHTED MEAN DIAMETER 
C DEL32 IS INCREMENT FOR D32 

C DIST IS SURFACE WEIGHTED LOG-NORMAL SIZE DISTRIBUTION 

C M IS THE MAXIMUM NUMBER OF POINTS IN Q-AREAY FOR MONODISPERSE 
C 

PRINT 240 

240 FORMAT(3X,'D32',SX,'DMIN',TX,'DMAX',5X,'NMBR',6X,'DlOX',7X,'D32X' 
1,7X,'QBAR'/) 

C 

D32=SD32 

DEL32*SDL32 

C NUMBER OF POINTS FOR PLOTS 2 AND 4 IS NPTS2 
NPTS2=ND-H 
DO 17 NN=1,ND 
DMIN=D32/SIGMA**3. 

DMAX=D32»SIGMA»*2.5 
AS=ALOG(SIGMA1**2. 

DG=EXP(ALOG(D32)-2.5*AS) 

QSUM=0.0 
DSUM=0.0 
XSUM=0.0 
DHOLD=0.0 
DNSUM=0.0 
XNSUM=0.0 
NMBR=0 
MM=M-2 

DO 13 N=1,MM 

IF {DQ(N,1).LT.DMIN) GO TO 13 
IF (DQ(N,1).GT.DHAX) GO TO 15 
IF {NMBR.GT.O) GO TO 242 
DMINZ=DQ(N,1) 

242 CONTINUE 

DD=DQ(N+1,1)-DQ(N,1) 

AA=(ALOG(DQ(N+l,l)/DG))**2/(2.0*AS) 
DIST=EXP(-AA)*DQ(N+1,1)/SQRT(6.2832*AS) 
QSUM=QSUM+(DIST»DQ(N+l,2)+DHOLD*DQ(N,2))*DD/2.0 
DSUH=DSUH-^(DIST*DQ(N+l,l)+DHOLD*DQ(N,l) )»DD/2.0 
XSUM=XSUM*(DIST+DKOLD)*DD/2.0 

XNSUM=XNSUM+-(DIST/DQ(N-^1,1)**2+DHOLD/DQ(N,1)*»2)/2.0*DD 
DNSUM-DNSUM-t-(DIST/DQ (N-^l, 1)+DKOLD/DQ (N, 1) )/2.0»DD 
XHOLD-D0(N-H,l) 

DHOLD=DIST 
NMBR-MMBR-t-1 
13 CONTINUE 
15 CONTINUE 

QBAR=OSUM/XSUM 
D32X=DSUM/XSUM 
D10X=DNSUM/XNSUM 
ALLQ(NN+1,2,NI)=D32X 
AliQ(NN+l,4,NI)=QBAR 
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DIAH32(NN+11=D32 

PRINT 241,D32,DMINZ,XHOLD,NMBR,D10X,D32X,QBAR 
241 F0RMAT(1X,3(F6.3,5X),I3,5X,2(F6.3.5X),P7.4,5X) 
D32=D32*DEL32 
17 CONTINUE 
RETURN 
END 
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APPENDIX D. MACRO FOR LOTUS 1-2-3 RELEASE 4 


On the following pages is a macro for Lotus* 1-2-3* Release 4 for Windows^M. 
When the output data file from ’'miescat3.exe" was imported into a 1-2-3 spreadsheet all 
the data would enter as a single column due to text blocks between sets of data. 

Therefore, data parsing was necessary to split the data and enter it into the appropriate 
number of columns. Since this operation needed to be performed every time the 
’'miescat3.exe" data was to be used it was advantageous to create a macro to automate this 
task. The macro was recorded as the data parsing was done manually for one entire file, 
then the macro was edited until it was streamlined and all the bugs were fixed. It is 
included here specifically for those who must also work with the ’'miescai3.exe'' output. 

It has been subsequently discovered that by using the File Open, [miescat3.exe 
output filename]. Combine, Formatted Text combination that the data will be imported 
properly. However, all of the text blocks including data labels will be missing. 

Therefore, the use of this macro may be the best option. 
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APPENDIX E. MIE CODE "MIESCAT6.FOR” 


To simplify the data reduction process (i.e,, finding the best combination of 
complex index of refraction m, geometric standard deviation o,. and Sauter mean 
diameter measured transmittances) the Mie Code in Appendix B was further 

modified. Each modification resulted in the program getting assigned a new number such 
as "miescat4.for''. Since the input data file and the compiled executable version usually 
also changed they were both given the new number. This particular program allows 
iteration starting values, step sizes, and the number of loops to be input in the data file for 
m and o,. The program outputs the current values of m and and the best fit rfjj as well 
as the Sum of the Squared Deviations (SSD) value for the fit between the InT ratios and 
the Q ratios. 
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BV KENNETH CASHDOLLAK, 1976 
REVISED BV THOMAS WELDON, JULY, 1977 
BASED ON PROGRAM HICOEF BV C D LITTON AND Z J FINK 
REVISED BV LT JOHN MANSER. USN 1994-1995 (PLOT 
FUNCTIONS REMOVED, UPDATED FOR MS-FORTRAH 4.1, 
DATA BLOCK WITH WAVES i T PUT INTO EXTERNAL FILE, 
DATA REDUCTION ITERATIONS ADDED) 


"MIESCAT6.EXE'' CALCULATES EXTINCTION COEFFICIENTS AND RATIOS FOR 
LIGHT TRANSMISSION THROUGH DUST PARTICLES WITH COMPLEX REFRACTIVE 
INDEX AND A LOG-NORMAL SIZE DISTRIBUTION AND OUTPUTS THE SUM OF 
THE SQUARED DEVIATIONS BETWEEN THE QBAR AND LN(T) RATIOS 


«*«>*«***«*** INPUT DATA FILE "HIESCAT6.DAT" FORMAT ••••»•••**•♦» 


VARIABLE NAMES AS IN PROGRAM 


C ROW 1....T1,T2,T3 

C ROW 2....WAV(1),WAV(2),WAV(3) 

C ROW 3....NUMR.STEPR 

C ROW 4....NUMS.STEPS 

C ROW 5....RDUST.RMED 

C ROW 6....L2,L3,DS,DEL0S 

C ROW 7_ND,SIGMA,SD32,SDL32 


FORMATS AND VARIABLE DESCRIPTION 

ROW ONE_F6.4(TRANSMITTANCE FOR WAVELENGTH 1) , 2X, F6.4(TRANSH 

ITTANCE FOR WAVELENGTH 2),2X,F6.4(TRANSMITTANCE FOR 
WAVELENGTH 3) 

ROW TWO_F6.4(WAVELENGTH 1),2X,F6.4(WAVELENGTH 2),2X,F6.4 

(WAVELENGTH 3) (ALL IN MICRONS) 

ROW THREE..I3(NUMBER OF ITERATIONS ON ROUST),2X,F4.3(STEP SIZE 
OF ITERATIONS). 

ROW FOUR...13(NUMBER OF ITERATIONS ON SIGMA),2X,F4.3(STEP SIZE 
OF ITERATIONS). 

ROW FIVE...2F6.3(COMPLEX REFRACTIVE INDEX OF DUST PARTICLES), 
2X,F6.3(REFRACTIVE INDEX OF THE MEDIUM). 

ROW SIX....2I3(L2,THE NUMBER OF INCREMENTS TAKEN AT A SET DELTA) 
(L3, THE NUMBER OF TIMES DELTA IS DOUBLED),2F6.3 
(STAR DIAMETER AND INITIAL VALUE OF DELTA). 

ROW SEVEN..I3(NUHBER OF INCREMENTS TAKEN FOR LOG-NORM),F4.2 
(SIGMA),2F5.2(STARTING DIAMETER AND DELTA). 



COMPLEX R,RDUST,RDUSTl,DI,DELDI,RX,RIHDX,DIS,DELIS 
COMPLEX COEFR,DPR(1000),XJR(1000),DJR(1000),PR(1000) 
DIMENSION XY(IOOO),WAV(3) 

DIMENSION DP(IOOO),XK(1000),XJ(1000),DJ(1000),P(1000) 
DIMENSION RATl(lOO),RAT2(100),RAT3(100) 

DIMENSION RAT4(100).RATS(100),RAT6(100),SSD(100) 
COMPLEX A(IOOO),B(1000),C(1000) 
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COMPLEX ALPH,XH(1000).DH(1000),E(1000),DE(1000) 
COMMON/ADQ/DQ(500,2) 

COMMON ALLQ(100,4,3) 

COMMON/KPLT/RINDX(3),WAVL(3) 

COMMON DIAM32(100) 


OPEN (UMIT=1, FILE='MIESCAT6.DAT') 

READ (1,103)T1,T2,T3 

READ (1,103)WAV(1),WAV(2),WAV(3) 

103 FORMAT(F6.4,2X.F6.4,2X,F6.4) 

READ (1,104)HUMR,STEPR 
READ (1,104)NUMS,STEPS 

104 FORMAT(I3,2X,F4.3) 

C DO 8 NDATA=1,10 

C READ (1,105)NDAT 

C 105 FORMAT(I3) 

C IF(NDAT.E0.999) GO TO 9 

NDAT=001 
C 

C RDUST IS COMPLEX REFRACTIVE INDEX OF DUST PARTICLES, WITH MINUS 
C RHED IS REAL REFRACTIVE INDEX OF MEDIUM 
READ (1,32)RDUST,RMED 
32 FORHAT(2F6.3,2X,F6.3) 

READ (1,34)L2,L3,DS,DELOS 
34 FORMAT(2I3,2F6.3i 

READ (1,12)ND,SIGMA,SD32,SDL32 
12 FORMAT(I3,F5.2,2F5.2) 

RDUST1=RDUST 
PRINT IS 

18 FORMAT!' EXP. TRANSMITTANCES:') 

PRINT 20,WAV(1),T1,WAV(2),T2,WAV(3),T3 

20 F0RMAT(3(' T(',F6.4,',F6.4)/) 

PRINT 21 

21 FORMAT!' BEST LOG-NORMAL FIT VALUES:') 

PRINT 22 

22 FORMAT!8X,'M',11X,'SIGMAG',7X,'D32',10X,'SSD') 

DO 9 INDEX1=1,NUMS 

DO 8 IHDEX2=1,NUMB 

C RX IS RELATIVE REFRACTIVE INDEX, (N-IK)/H 

RX=ROUST/RHED 

C IMAGINARY PART OF R CHANGED TO POSITIVE TO CONFORM TO EXPANSION 
C FORMULA IN PROGRAM 

R=CONJG!RX) 

RINDX(1)=RX 

RINDX(2)=RX 

RINDX!3)=RX 

C 

C L2*L3 = NUMBER OF SIZE INCREMENTS 

C INITIAL DIAMETER IS DS, INITIAL VALUE OF DELTA IS DELDS 


C INITIAL VALUE OF COMPLEX DIAMETER IS DIS 
DIS=CMPLX(0.0,DS) 

DELIS-CMPLX!0.0,DELDS) 

C INITIAL VALUES FOR ARRAY 
DO 51 K=l,3 
DO 52 J=l,2 
ALLQ!1,J,K)=0.0 
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ALLQ(1,J+2,K)=0.0 
DO 53 1= 2,100 
ALLQ(I,J,K)=100.0 
ALLQ(I,J+2,K)-1.0 

53 CONTINUE 
52 CONTINUE 
51 CONTINUE 

DO 54 1=1,100 

DIAM32(I)=0.0 

RAT1(I)=0.0 

RAT2(I)=0.0 

RAT3(Il-0.0 

54 CONTINUE 

DO 4 NI=1,3 
C 

C PRINT 106,NDATA 

C 106 F0RMAT(/' DATA SET NUMBER',13) 

C PRINT 35 

C 35 FORMAT(' HIE COEFFICIENTS FOR MONODISPERSE DUST PARTICLES') 

C PRINT 41,RDUST,RHED,RX 

C 41 F0RMAT(' RDUST=',2F7.4,/' RMED=',F7.4,/' REL. REF. INDEX=',2F7.4) 
C PRINT 104,B 

C 104 FORMAK' R=',2F7.4) 

WAVE=WAV(NI) 

C PRINT 30,WAVE 

C 30 format;' WAVELENGTH =', F7.4,2X,'MICRONS'/) 

C PRINT 25 

C 25 FORMAT(6X,'A',12X,'D',13X,'QABS',11X,'QEXT'/) 
WAVL(N1)=WAVE*10000.0 
C 
C 

C INITIAL VALUES SET 

C MAX VALUE FOR D IS (1000*WAVE) / (1.3*PI*RMED) = 244*WAVE/RMED 
C D IS REAL DIAMETER IN MICRONS,DELD IS INCREMENT 
C DI=I»D IS COMPLEX DIAMETER. DELDI IS INCREMENT 

C X IS REAL SIZE PARAMETER X= PI*DIAMETEB»RHED/WAVELENGTH 
C ALPH IS COMPLEX SIZE PARAMETER, ALPH=I*X 

D=DS 

DELD=DELDS 

DI=DIS 

DELDI=DELIS 


DO 3 J=1,L3 

DO 2 1=1, L2 

X=3.14159»D/WAVE*RMED 

ALPH=3.14159*DI/WAVE*RMED 

XJ(1)= SIN(X)/X 

XJ(2)= SIN(X)/(X**2)- COS(X)/X 

XJR{1)=CSIN{R*X)/(R*X) 

XJRC2)=CSIN(R*X)/( CR*X)->*2)-CCOS(R*X)/(R»X) 
XH(1)“CEXP(ALPH)/ALPH 

XH{2)=((-CEXP(ALPH))/X)-((0.,1.)*CEXP(ALPH)/(X**2)) 
DJ(1)= COS(X)/X- SIN(X)/(X»*2) 
DJR(l)=CCOS(R*X)/(B*X)-CSIN(R*X)/{(R*X)»*2) 
DH(1)=CEXP( ALPH)/X-t-((0. , 1. ) *CEXP(ALPH)/( X*»2) ) 
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QEXT=0.0 

QSCA=0.0 


C LI = ITERATIONS FOR SPHERICAL FUNCTIONS, L1=1.I*X 
L1=1..3*X 

IF(Ll.LE.lO) Ll=10 
DO 1 L=1.L1 
CO-(FLOAT(L) )/(2*L+l.) 

COE=(FLOAT(L+1))/(2*L+1) 

COEF=(FLOAT(2*L+1))/X 
COEFR=(FLOAT(2*L+1))/(R*X) 

XJ(L+2)=COEF*XJ(L+l)-XJ(L) 

XJR(L+2)=C0EFR*XJR(L+1)-XJR(L) 

XH(L*2)-COEF*XH(L+l)-XH(L) 

DJ(L-*-l) =CO»XJ( L) -COE*XJ (L+-2 ) 

D JR (L-fl) =CO* XJR (L)-COE* XJR (L-t-2) 

DH(L*l)=CO*XH(L)-COE*XH(L+2) 

P(L)=X*XJ(L+1.) 

PR(L)=(R*X)*XJR(L+1) 

DP(L)=XJ(L+-l)-t-X*DJ(L+l) 

DPR(L)=XJR<L-t-l) + (R»X)*DJR{L-t-l) 

E(L)-X*XH(L+1) 

DE(L)=XH(L+1)+X*DH(L+1) 

A(L)=((DPR(L)*P(L))-(R*PR{L)*DP(L)))/((DPR(L)*E(L))-(R*PR(L)*DE(L 

iin 

B(L)=((R*DPR(L)*P(L))-(PR(L)*DP(L)))/((R*DPR(L)*E(L))-(PR(L)*DE(L 

C(L)-A(L)+BCL) 

GAMMA=CABS(A(L)) 

BETA=CABS(B(L)) 

XK1L) = (FL0AT(2*L+1))*(2./(X**2))*{ (GAMHA**2)-t-(BETA**2)) 
XY(L)=(FL0AT(2*L+:))*(2./(X**2))*(REAL(C(L))) 

QEXT=QEXT-*-XY{L) 

QSCA=(3SCA*XK(L) 

1 CONTINUE 

QABS=QEXT-QSCA 

DQ(M,1)=D 

DQ(H,2)=QEXT 

C PRINT 22,X,DQ(H,1),QABS,DQ(H,2) 

C 22 FORMAT(1X,F0.3,6X,F7.3,9X,F6.3,9X,F6.3) 

C 

C D AND QEXT FOE HONODISPERSE PUT INTO PLOTTING ARRAY 
C EVERY FOURTH POINT CALCULATED IS POT INTO ARRAY- ALLQ 


XM»H 



XHINT=MINT 

IF (AREL.NE.XHINT) GO TO 60 
IF (D.GT.1.20) GO TO 60 

C NUMBER OF POINTS FOR PLOTS 1 AND 3 IS NPTSl 

NPTS1=HINT 
AIXQ(MINT+1,1,NI)=D 
ALLQ(HINT-t-l, 3 ,NI)=QEXT 
ALLQ(HINT+1,1,NI)=D 
ALLQ (MINT+-1,3 , NI ) =QEXT 
60 CONTINUE 

C INCREMENT SIZE PARAMETERS D=DIAM£TER 

D=D-t-DELD 
DI=DI+DELDI 
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M=M+1 

2 CONTINUE 

C PRINT 36 

36 FORMAT (IK ) 

DELD=2.0*DELD 

DELDI=2.0»DELDr 

3 CONTINUE 

C PRINT 106.NDATA 

C PRINT 235 

C 235 FORMAT(' MIE COEFFICIENTS FOR LOG-NORMAL DUST DISTRIBUTION') 

C PRINT 236,SIGMA 

236 FORMAT (' SIGMA-',F4.2) 

C PRINT 237,RX 

237 FORMATC RELATIVE REFRACTIVE INDEX-' ,2F7.4 ) 

C PRINT 30,WAVE 

CALL QAVG(M,HI,SIGMA,NPTS2,SD32,SDL32,ND) 

4 CONTINUE 

C PRINT 106,NDATA 

7 CONTINUE 

C PRINT 500 

C 500 FORHAT[' QBAR RATIOS FOR 3 WAVELENGTHS') 

C PRINT 236,SIGMA 

C PRINT 237,RX 

C PRINT 509,WAV(1),T1,WAV(2),T2,WAV(3),T3 

C 509 F0RMAT(' EXP. TRANSHITTANCES:',3(5X,'T(',F6.4,')=',P6.4)) 

C PRINT 510,WAV(3),WAV(1),WAV(3),WAV(2),WAV(2),WAV(1) 

C 510 FORHAT(/2X,' D32 ',3(5X,'QBAR(',F6.4,')/QBAR(',F6.4,') ' ),5X,'SUM 
C +0F SQUARED DEVIATIONS'/) 

K-ND+1 

DO 600 N=2,K 

RAT1(N)-ALLQ(N,4,3)/ALLQ(N,4,1) 

RAT2(N)=ALLQ(N,4,3)/ALLQ(N,4,2) 

RAT3(N)=ALLQ(N,4,2)/ALLQ(N,4,1) 

RAT4(N)=(RAT1(N)-LOG(T3)/LOG(T1))»*2 
RAT5(N)=(RAT2(N)-LOG(T3)/LOG{T2))»*2 
RAT6(N)=(RAT3(N)-LOGCT2)/LOG(T1))»»2 
SSD(N)-RAT4(N)+RAT5(N)+PAT6(N) 

C PRINT 650,DIAM32(N),RAT1(N),RAT2(N),RAT3(N),SSD(N) 

600 CONTINUE 

C 650 FORMAT(lX,F6.3,13X,F7.4,23X,F7.4,23X,F7.4,23X,F9.4) 

C PRINT 350,NDAT 

C 350 FORMAT (/'•*•••••*•••••••••••••«*••» END OF DATA SET NUMBER',13/) 

SMALLEST-2 
DO 700 N=2,K 

IF (SSD(N).LT.SSD(SMALLEST)) SMALLEST-N 

700 CONTINUE 

PRINT 701,RDUST,SIGMA,DIAK32 (SMALLEST) ,SSD(SMALLEST) 

701 FORMAT (IX,2F7.4,'I',5X,F4.2,6X,F5.2,6X,F13.8) 

NDAT-NDAT+1 

RDUST-RDUST+STEPR 

8 CONTINUE 
ROUST-RDUSTl 
SIGMA-SIGMA+STEPS 

9 CONTINUE 
END 


SUBROUTINE QAVG(M,NI,SIGMA,NPTS2,SD32,SDL32,ND) 


B'i KENNETH CASHDOLLAR, 1976 





C QAVG CALCULATES D32 AND QBAB FOR A GIVEN QEXT VS D AND 
C A GIVEN LOG-NORMAL SIZE DISTRIBUTION 


C INPUT DATA = ND,SIGMA,D32,DEL32 


COMMON ALLQ(100,4,3) 

COMMON/ADQ/DQ(500,2) 

COMMON DIAM32(100) 

C DQ(I,1)“D, DQ(I,2)=QEXT 

C NI=WAVELENGTH OR REFRACTIVE INDEX CHANGE 

C ND IS NUMBER OF D32,(3BAR CALCULATIONS 
C SIGMA IS THE GEOMETRIC STANDARD DEVIATION 
C D32 IS THE SURFACE WEIGHTED MEAN DIAMETER 

C DEL32 IS INCREMENT FOR D32 

C DIST IS SURFACE WEIGHTED LOG-NORMAL SIZE DISTRIBUTION 

C H IS THE MAXIMUM NUMBER OF POINTS IN Q-ARRAY FOR MONODISPERSE 

C PRINT 240 

C 240 FORMAT(3X,'D32',8X.'DMIN',7X,'DHAX',5X,'NMBR',6X,'DlOX',7X,'D32X' 
C 1,7X,'QBAR'/) 

C 

D32=SD32 

DEL32=SDL32 

C 

C NUMBER OF POINTS FOR PLOTS 2 AND 4 IS NPTS2 
NPTS2=ND+1 
DO 17 NN=1,ND 
DHIN=D32/SIGMA»»3. 

DKAX=D32»SIGMA**2.5 
AS=ALOG(SIGMA)* * 2. 

DG=EXP(ALOG(D32)-2.5»AS) 

QSUM=0.0 
DSUM-0.0 
XSUM=0.0 
DHOLD=0.0 
DNSUM=0.0 
XNSUM=0.0 
NMBR=0 
KM=M-2 

DO 13 N=1,MM 

IF (DQ(N,1).LT.DMIN) GO TO 13 
IF (DQ(N,1).GT.DMAX) GO TO 15 
IF (NMBR.GT.O) GO TO 242 
DMINZ-DQ(N,1) 

242 OONTINUE 

DD=DQ(N+1,1)-DQ(N,1) 

AA=(AL0G(DQ(N-h1,1)/DG) )**2/(2.0*AS) 
DIST-EXP(-AA)»DQ(N*1,1)/SQRT(6.2832*AS) 

0SUM=QSUM+(DIST*DQ(N+1,2)+DH0LD*DQ(N,2))*DD/2.0 
DSUM=DSUM+(DIST*DQ(N-H,1)+DHOLD*DQ(N,1))*DD/2.0 
XSUM=XSUM+(DIST+DHOLO)»DD/2.0 

XNSUM=XNSUM+(DIST/DQ(N*1,1)»»2+DHOLD/DQ(N,1)**2)/2.0*DD 
DNSUM“DNSUM*(D1ST/DQ(N*1,1)+DHOLD/DQ(N,l))/2.0*DD 
XH0LD=DQ(N+1,1) 

DHOLD-DIST 
HMBR=NHBR+1 
13 CONTINUE 
15 CONTINUE 
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QBAS=QSUH/XSO>l 
D32X=DSUM/XSUH 
DIOX-DNSUM/XNSUH 
ALLQ(NN+1,2,NI)”D32X 
ALLQ(NN*1,4,NI)=QBAR 
DIAH32(NN+11-D32 

PRINT 241,D32,DMINZ,XHOLD,NMBR,D10X,D32X,«BAR 
241 rORMAT(lX,3(F6.3,5X),I3,5X,2(F6.3,5X),F7.4,5X) 
D32=D32+DEL32 
17 CONTINUE 
RETURN 
END 
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APPENDIX F. LINEAR INTERPOLATION IN LOTUS 1-2.3 RELEASE 4 


A formula was written to perform linear interpolation in Lotus® 1-2-3® Release 4 
for Windows™. The formula used the diameter value in column B (cell address BIO 
below) to choose the correct extinction coefficient Q from a separate 1-2-3 file named 
"qext.wkd". The file "qext.wk4" had a named range QEXT containing seven colunms of 
data: diameter and six columns of Q's for different wavelengths and indices of refraction. 
The formula below contains the offset for the first column of Q’s, 

@INDEX(«qext-wk4»$QEXT.l,@MATCH(B10.«qext-wk4»$QEXT,!))+(( 
@INDEX(«qext.wk4»$QEXT,l,@MATCH(BiO.«qext.wk4»$QEXT.l)+l)- 
@INDEX(«qext-wk4»$QEXT,l,@MATCH(B10.«qext-wk4»$QEXT,l)))/(@ 
rNDEX(«qext.wk4»$QEXT.O.@MATCH(B10,«qext.wk4»$QEXT,l)-(-l)-@I 
NDEX(«qext.wk4»$QEXT,0,@MATCH(B10,«qext.wk4»$QEXT.l)))*(B10 
-@rNDEX(«qext,wk4»$QEXT,0.@MATCH(B10,«qext.wk4»$QEXT,l)))) 
Linear interpolation can best be described by an example. Suppose d, and are 
increasing diameters with corresponding extinction coefficients Q, and Q^. Then the Q 
for a diameter d between d, and is: 


The reader is directed to the User’s Guide for further information and the ©function help 
in the software for an explanation of ©INDEX and ©MATCH, 
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Figure 11. PDPA Data from PS-1 Firing 3/2/95. 
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Diameter Distribution 


Bin Diaieter Count Bin Diaieter Count Bin Diaieter Count Bin Dianeter Count Bin Dianeter Cm 


0 0.500 1 8 1.320 2 16 2.H0 0 2< 2.960 0 32 3.780 

1 0,603 1 9 l.t23 0 17 2.243 0 25 3.063 0 33 3.883 

2 0.705 0 10 1.525 0 18 2.345 0 26 3.165 2 34 3.985 

3 0.808 0 11 1.628 3 19 2.448 0 27 3.268 0 35 4.088 

4 0.910 1 12 1.730 1 20 2.550 2 28 3.370 1 36 4.190 

5 1.013 0 13 1.833 0 21 2.653 1 29 3.473 1 37 4.293 

6 1.115 0 14 1.935 1 22 2.755 1 30 3.575 1 38 4.395 

7 1,218 0 15 2.038 0 23 2,858 2 31 3,678 0 39 4.498 


Bin Diaieter Count Bin Diaieter Count Bin Diaieter Count Bin Diaieter Count Bin Diaieter Co 


40 4.600 1 62 6.856 0 84 

41 4.703 0 63 6.958 0 85 

42 4.805 0 64 7.061 0 86 

43 4.908 0 65 7.163 0 87 

44 5.010 0 66 7.266 0 88 

45 5.113 0 67 7.368 0 89 

46 5.215 2 68 7,471 1 90 

47 5.318 0 69 7.573 0 91 

48 5.421 1 70 7.676 0 92 

49 5.523 0 71 7.778 0 93 

50 5.626 1 72 7.881 0 94 

51 5.728 1 73 7.983 0 95 

52 5.831 1 74 8.086 , 0 96 

53 5.933 0 75 8.188 0 97 

54 6.036 0 76 8.291 0 98 

55 6.138 0 77 8.393 0 99 

56 6.241 0 78 8.496 0 100 

57 6.343 0 79 8.598 0 101 

58 6.446 0 80 8.701 0 102 

59 6.548 0 81 8.803 0 103 

60 6.651 0 82 8.906 0 104 

61 6.753 0 83 9.008 1 105 


ui ui ui 

9.111 1 106 11.366 0 128 13.621 

9.213 0 107 11,469 0 129 13.724 

9.316 0 108 11.571 0 130 13.826 

9.418 0 109 11.674 1 131 13.929 

9.521 0 110 11.776 0 132 14.031 

9.623 0 111 11,879 0 133 14.134 

9.726 0 112 11,981 0 134 14.236 

9.828 0 113 12.084 0 135 14.339 

9.931 0 114 12.186 0 136 14.441 

10.033 0 115 12.289 0 137 14.544 

10.136 0 116 12.391 0 138 14.646 

10.238 0 117 12,494 0 139 14.749 

10.341 0 118 12.596 0 140 14.851 

10.444 0 119 12.699 0 141 14.954 

10.546 0 120 12.801 0 142 15.056 

10.649 0 121 12.904 0 143 15.159 

10.751 0 122 13,006 0 144 15.262 

10.854 0 123 13.109 0 145 15.364 

10.956 0 124 13.211 0 146 15.467 

11.059 0 125 13.314 0 147 15.569 

11.161 0 126 13.416 0 148 15.672 

11.264 0 127 13.519 0 149 15.774 


Figure 11 Confd, PDPA Data from PS-1 Firing 3/2/95. 
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Due to a computer hard drive problem, this data was incapable of being stored or 
further analyzed. It was interpreted from the distribution plot and the mean diameters and 
recorded by hand. PDPA bin diameter guesses that resulted in the calculated arithmetic 
mean and Sauter mean rfjj were then used and are presented below. 


PDPA bin diameter Number 

Arithmetic Mean (DIO) 

1= 2.3 urn 

1.013 nm 2 

Sauter Mean (D32) 

= 3.7 p.m 

3.063 |im 1 

PVC not calculated 


4.293 (am 1 

Attempts 

851 


Valid 

4 


% Valid 

0 


Figurel2. PDPA Data from PS-1 Firing 3/3/95. 




Figure 13. PDPA Data from PS-1 Firing 3/4/95. 


























PS-1 fIRING 2/25/95 



TIME in seconds 


Figure 14. LABTECH Notebook Data Acquisition Program Plot of PS-1 Firing 2/25/95. 
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0-2 sec 
8-78-8.82 sec 
no light 


Diode Outputs (v) 

0.850Q 0.5461 0.3130 

2.5443 4,0226 2.1436 

0.8406 0.8894 0.9525 

0.6600 0.7131 0.8822 

0.0958 0.0533 0.0557 


Figure 15. Normal Plot ofPS-1 Firing 2/25/95. 
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0-2 sec 
5.4-6.9 sec 
no light 


Diode Outputs (v) 

0.8500 0.5461 0.3130 


8.3110 

6-1475 

0.6270 

0.7184 


7.7200 

6.2676 

0.6700 

0.7940 


6.3999 

9.9938 

0.8060 

1.6425 


Figure 16. Bare Light (No Probe) PS-1 Firing 3/5/95. 
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Number & Mass P(D) & F{D) vs. D (microns) 
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Figure 18. MATLAB Distribution Plots for "Small Particles". 
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Figure 19. Scanning Electron Microscope Pictures of AljOj Collected in the Probe. 
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